
m the mferestof early m^mSf dfs- 
ssminetWi of ^^sources Survey 
"Cijram ^formation an'^rthout i»£»jjtv 
J*se mad# thir«of. ir 




***«»»** 


E7. 4- 1 0-7 !•• 
c^- njtrc 


KANSAS ENVIRONMENTAL AND 
RESOURCE STUDY: A GREAT PLAINS 
MODEL 

MAY 1974 


Type III Final Report for the 
Period August 1, 1972 through May 15, 1974 

-• t 


Original photography may. be. giirchased IIQIIU 
EROS Data Center 

10th and Dakota Avenue' , 

Sioux Falls, SO. &7JM 

*«4^A>£jyt 



Prepared for: 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 

Contract No. NAS 5“21822, Task 6 



THE UNIVERSITY OF KANSAS CENTER FOR RESEARCH, INC. 

2385 Irving Hill Rd— Campus West Lawrence, Kansas 66044 /F* r^r f! ^ ;<i J '».;•)•) 

AiinjfS 1 s: - i 8i77i 
' ^ rr Y} ^ 


KANSAS ENVIRONMENTAL AND RESOURCE STUDY: 
A GREAT PLAINS MODEL 


Ground Pattern Analysis in the Great Plains 


F. T. Ulaby, Principal Investigator (Acting) 

University of Kansas Center for Research, Inc. 

Remote Sensing Laboratory 
c/o Space Technology Center 
Nichols Hall 

2291 Irving Hill Drive“Campus West 
Lawrence, Kansas 66045 

(E74- 1 Q7 1 4) GROUND PATTERN ANALYSIS IN N74-31795^ 

THE GREAT PLAINS Final Report, 1 Aug. 

1972 - 15 May 1974 (Kansas Univ. Center 

■for Research, Inc.) 84 p HC $7.25 Unclas 

' -n r - -- CSCL 08B G3/13 00714 J 

May 1974 

Type III Final Report for the period August 1, 1972 through May 15, 1974 
Report No. 2266—1 1 


O^ral photography may. be. fiUBtesi IBW 

enos Data Center 
Win and Dakota ton* 

Sioux Falls, SD. ^ 


Prepared for: 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
GODDARD SPACE FLIGHT CENTER 
GREENBELT, MARYLAND 20771 

Contract No. NAS 5— 21 822, Task 6 



Final ERTS-A User 
Investigation Report 


Type III Final Report 

for the Period August 1, 1972 through May 15, 1974 


NASA Contract NAS 5- 21822 


Title of Investigation: Ground Pattern Analysis in the Great Plains 
ERTS-A Proposal No.: 60“8 
Task Number: 6 

Co-principle Investigators: John C. Davis and Fawwaz T. Ulaby 
NASA-GSFC PI ID No.: UN 657 


Report Prepared by: 


Report Approved by: 




James L. McNaughton 
Research Engineer 





Fawwaz T. Ulaby 
Co-principle Investigator 



TECHNICAL REPORT STANDARD TITLE PACE 


j. kvcipivnl't Catolog Ho. 

5”. Report t>at* 

May 1974 


6 , P*fformrng Or 90 ** *at<or> Cod* 


8. Pe r for ni no Orgom x otton Report No- 
2266-11 


10. Work Unit No. 


H- Contract Of Grant No. 

NAS 5-21822 

13. Type of Report and Period Covered 

Type III Final Report for 
the period Aug, 1, 1972 
to May 15, 19 74 

,14. Sponsoring Agency Code 


. Report No. 2. Government Accession No. 

2266.-1 1 - 1 

T. T»tie «nj SwLr.tl* |< ansas Environmental and Resource 

Study: A Great Plains Model 

Ground Pattern Analysis in the Great Plains 


7. Author(t) 

James L. McNaughton and Fawwaz T. Ulaby 


9. Performing Orgoni jaf ion Name ortd Address 

University of Kansas Center for Research, Inc. 
c/o Space Technology Center 
2291 Irving Hill Dr. -Campus West 
Lawrence, Kansas 66045 _ 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 
G. R. Stonesifer, Code 430 . 

15. Supplementary Notes 


None 

16. Abstract 

The purpose of this investigation was to locate large-scale geologic ground 
patterns and anomalies in Kansas using Earth Resources Technology Satellite 
(ERTS-1) imagery and optical data processing techniques. ERTS-1 sample areas 
were chosen and used as the input to an optical data processing system. The op- 
tical processor forms the diffraction pattern of the sample area and quantitative 
information obtained from the diffraction pattern is used to charcterize each 
sample area. This characterization of the ERTS sample areas is shown to work 
very well in discriminating between ground pattern categories in Kansas. The 
method of analysis used here is general and could be used to process various 
types of images. 


17. Key Words (S. lccted by Aofiior(i)) 

Optical data processing; Diffraction 
pattern sampling; Spatial frequency 
analysis; Large-scale geologic ground 
pattern identification 


18. Dittribution Statement 

Original, photography may be Rurcbasfid IiQflU 
EROS Oata Center 
10 th and Dakota Avenue 
Sioux Fails* SD 57138 


19. Security Ctosuf- (of tfiit report) 

20. Security Clatsif. (of *!>»» page) 

21 . No. of Paget 

Unclassified 

Unclassified 

75 



1 1 







PREFACE 


The objectives of this investigation are: 

A. Mapping the surficial geology of selected sites in Kansas from 
multispectral imagery, and identification of anomalous patterns; 

B. Search for large-scale ground patterns by spatial frequency analysis. 

An optical processing system is used in this investigation to produce ground 
pattern spatial frequency and orientation information using ERTS-1 imagery as input, 
interpretation of the information is done with respect to known geologic features and 
other cultural and vegetation features. Appropriate data processing schemes are used 
to derive numerical descriptors of these features and these descriptors, in turn, are 
used in pattern recognition schemes. This investigation will then provide a mapping 
technique of large-scale geologic ground patterns as well as other large-scale ground 
patterns in Kansas. 

The manual interpretation of the spatial frequency and orientational informa- 
tion derived from the ERTS-1 imagery demonstrates that this information characterizes 
the samples areas very well. The results described here, based on this characterization 
of the sample areas, show that optical diffraction analysis can be used to accurately 
discriminate between large-scale ground patterns in Kansas. 
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1,0 INTRODUCTION 


The purpose of this investigation was to determine if spatial frequency anal- 
ysis of Earth Resources Technology Satellite (ERTS-1) imagery could provide an un- 
biased means of estimating the character and nature of large-scale ground patterns 
in the image. Specifically, this investigation involved the determination of spatial 
frequency and orientational information from large-scale ground patterns in Kansas 
using ERTS-1 imagery as the input to an optical data processing system. Since the 
physiography and geology of Kansas is well-known, this method would provide an 
objective means of identifying and classifying ground patterns and locating anomalous 
patterns. 

The optical processing system used in this investigation forms the diffraction 
pattern of the input image transparency and provides data on the spatial spectrum 
of the image. It will be shown in this paper that this spatial frequency analysis of 
ERTS-1 images provides an accurate and unbiased means of automatically identifying 
large-scale ground patterns in Kansas. It is anticipated that the equipment, procedures, 
and methods described in this paper could be applied in general to obtain automatic 
pattern recognition using various types of images. 

This paper will describe the physiographic and geologic patterns that exist 
in Kansas; it will review the theory involved in diffraction pattern analysis; it will 
describe the optical processing system used in this investigation; it will review the 
spatial frequency and orientational information derived from this optical processor; 
it will state the accuracy with which these parameters predict the existence of large- 
scale ground patterns and anomalous patterns in Kansas; and finally, this paper will 
state some of the conclusions that can be made based on the results obtained during 
the course of this investigation. 

The procedure followed during this investigation is outlined below: 

(1) Determine the number and extent of various geologic/physiographic 
regions in Kansas from known geologic/physiographic information. 

(2) Select an adequate number of sample areas within each region. 

(3) Investigate the effect of snow cover on the ability of ERTS-1 imagery 
to display physiographic information. (To study this effect, equi- 
valent sample areas are chosen from imagery with and without snow 
cover.) 
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(4) Obtain data from the optical processor for the given sample areas. 

(5) Determine features of the optical processor data which will allow the 
categorization of the sample areas. 

(6) Investigate the ability of optical processor data to provide reliable 
physiographic classification using different classification algorithms. 

(7) Investigate the information content of various spatial frequency bands 
with respect to the investigation in (6) above. 

Eight regions of Kansas were chosen as most representative of certain dis- 
tinctive physiographic/geologic provinces. Ten sample areas were chosen from each 
of these regions. Each sample area is circular with a diameter of approximately 37 
kilometers. Optical processor data were taken twice for each of these sample areas 
once using an image when the area was snow-covered and once again using an image 
from the area when it was not snow-covered. 

Spatial frequency and orientational curves for 148 sample areas (80 non- 
snow and 68 snow samples) were obtained using the optical data processing system. 
Appropriate data processing schemes were used to derive numerical parameters from 
these curves, and these parameters are used as descriptors of the sample areas. These 
parameters were used to categorize the sample areas in terms of the different physio- 
graphic/geologic provinces previously chosen. 

It should be noted that the sample area sites were chosen where stream pat- 
terns provide a good expression of the geology or physiography of the geologic/physio - 
graphic province. Stream patterns, with respect to the scale of this experiment, are 
numerous, however, so this restriction did not severely limit the choice of sample 
sites . 

The images that were used in this experiment were obtained from NASA's 
ERTS-1 Multispectral Scanner (MSS) Subsystem. Each E RTS - 1 image covers an area 
of approximately lOOby 100 nautical miles. 

Categories and Sample Site Location 

Eight Physiographic regions in Kansas are shown in Figure 1 as adapted from 
Schoewe (1949). Each physiographic region is also a region of characteristic geology 
since each region contains similar landforms caused by similar geomorphic processes 
and each region also contains outcrops of the same dominant lithology and the same 
geologic age. An important note with respect to this investigation is that, in general, 
the land-use for each region is the same. Hence, each sample site was chosen in 
a well-defined physiographic/geo logic region of Kansas. 
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Ten sample areas were chosen In each of these categories and the sample 
area center points are also plotted in Figure 1 (Some sample areas are in or cross 
over into neighboring states but are still in the same geologic category). Each site 
has an alphanumeric designation. F - 1 refers to the first sample area in the Flint 
Hills; 0~6 refers to the sixth sample area in the Osage Plains; etc. When the image 
sample areas were analyzed with snow - cover they were assigned a number ten more 
than the non-snow number (0-6, no snow becomes 0-16, snow). The sample area 
in the High Plains region of Kansasdesignated "M n in Figure 1 was analyzed for 
all four MSS frequency bands (MSS 4, 5, 6, 7). This was done in order to determine 
the relative information content of each band. For all other sample areas, MSS 
band 5 (0.6-0. 7 microns) 70 mm positive transparencies were used in the analysis. 

It was judged that band 5 images contain the most information with respect to the 
objectives of this investigation. 

As mentioned previously, each sample area is approximately 37 kilometers 
in diameter and this diameter plus the size of the physiographic region controlled 
the spacing of the sample areas within each category. The sample areas were analyzed 
twice once with snow cover and once without snow cover. Most of the non-snow 
samples were from images recorded in the summer and fall of 1972. The images for 
samples with snow-cover were recorded in the winter of 1972/1973. Only 68 of the 
80 sample areas were analyzed with snow-cover, however, since cloud-free images 
were not available for 12 of the snow-covered sample areas. 

Detailed analysis of each sample area was performed and the results of this 
analysis were presented in Ulaby et al. (1973). For a discussion of the physiography 
and geology of Kansas see Appendix C. 
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Figure 1. Physiographic Regions of Kansas (Adapted from Schoewe , 1949 ) with Sample Site 
Location. 




2.0 BACKGROUND 


2.1 Theoretical Considerations 

This section is concerned with reviewing some of the theoretical background 
necessary to understand the basis of this investigation . The general area described 
here is referred to as optical data processing; the particular area with which this 
investigation is concerned is referred to as optical diffraction analysis or spatial 
frequency analysis. These terms arise because of the nature of the optical systems 
used. 

The basic element needed to describe the nature of an optical data processing 
system for optical diffraction analysis or spatial frequency analysis is a single positive 
lens and an input object illuminated by an incident plane wave. From basic physical 
optics we know that a positive lens illuminated by an incident plane wave converts 
this wave into a wave converging upon a point f, where f is the focal point of the 
lens. The amazing aspect of this simple lens system is Its ability to perform a two- 
dimensional Fourier transform. 

If a plane object is inserted into this simple system at a point d^ away from 
the lens as shown In Figure 2, the distribution in the focal plane of the lens Is pro- 
portional to the Fourier transform of the object's transmittance function (i.e. the com- 
plex function associated with the way in which the object transmits light). The 
illumination is assumed to be monochromatic and coherent. 

In particular, if 

t Q (x,y) is the complex amplitude transmittance 
function associated with the portion of 
the object illuminated in Figure 2, and 

UjCx^y,) is the complex amplitude distribution 
in the back focal plane of the lens, and 

A is the amplitude of the normally incident 
plane wave, 

and assuming that the portion of the object illuminated is smaller in extent than the 
clear aperture of the lens we have (Goodman (1968) ): 
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where 


h'/i) “ f 0 ( x /y) e x p “j 77 — (xx, +yy,) dxd y 0 ) 

-00 L 1 1 J 



Hence, the diffraction pattern in the back focal plane of the lens is proportional 
to the Fourier transform of the function which describes the object. 

Note that the terms x^/f\ and y-j/fX in Equation (1) are scale factors 
relating the scale of the frequency distribution in the transform plane (back focal 
plane) to the scale of the distribution in the object plane, (i.e. for any point 
(x^,yj) in the transform plane, the spatial frequency associated with this point may 
be calculated from the geometry of the system to obtain the spatial frequency com - 

*1 yi 

ponent values ( — — ) ). 

A mathematical description of the system used in this investigation is given 
in Appendix A. For a description of other optical processing configurations, the 
reader may consult Goodman (1968), Shulman (1970), or Preston (1972). 

The distribution of interest has been described in terms of the Fourier transform 
of the input amplitude transmittance function. Hence, well-known theories from 
Fourier analysis can be applied in relating the transform to the original input (image 
transparency in this case). 

Some of the aspectsof Fourier analysis as applied to this investigation were 
presented in Ulaby et al. (1973) and McCauley et al. (1974). The basics of Fourier 
analysis can also be found in a number of textbooks dealing with advanced mathe- 
matics as well as in a number of electrical engineering texts dealing with the spec- 
tral analysis of signals. 
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2.2 Other Investigations 

Basic physical optics and diffraction problems are discussed in Born and Wolf 
(1959), Goodman (1968), and Parrent and Thompson (1969). Optical systems similar 
to the one used in this investigation are discussed in Lee and Gossen (1971), and Read 
and Cannata (1974). Other current investigations which have used optical diffraction 
analysis for image analysis are: Stanley, Nienow and Lendaris (1969), and Lendaris 
and Stanley (1969), in which diffraction pattern sampling is used as a means to attain 
automatic target recognition; Davis and Preston (1972), McCuIlagh and Davis (1972), 
Pincus and Dobrin (1966), and Steckley (1972), in which optical diffraction analysis 
is used in geologic applications; and Read and Cannata (1974), in which one - dimen - 
sional data records are analyzed. Another investigation in which ERT5 imagery is 
analyzed by diffraction pattern analysis is Gramenopoulos (1973). Two investigations 
which illustrate the technique of Image processing by digital means are Haralick and 
Shanmugan (1973), and Andrews et al. (1972). 
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3.0 


EXPERIMENTAL ARRANGEMENT 


3.1 Spatial Frequency Analysis of Sample Areas 

Each sample area on the ERTS image is composed of ground patterns of 
characteristic spatial frequencies and orientations. The spatial frequencies here 
refer to variations of density on the image as a function of distance. Hence, if 
we decompose the image sample area into Its component frequencies and plot their 
intensity; and if we plot the 'strength' of preferred orientations in the image versus 
angle, we have described the sample area by a set of numbers. The purpose of the 
quantitative phase of this investigation was to determine how accurately this could 
be done • 

The optical data processing system used in this investigation accomplishes 
this decomposition of the sample area in the manner described in section 2.1 . A 
set of numbers is obtained which describe the spatial spectrum of the sample area 
and indicate the orientations of the large-scale ground patterns. Each of the sample 
areas on the ERTS images is used as the input to this optical processor. 

A block diagram of the optical processing system along with the pattern 
recognition portion of the system is shown in Figure 3. Figure 4 expands on the part 
of Figure 3 enclosed In dashed lines. These diagrams show the various steps involved 
in going from the ERTS sample area to a categorization of the sample area based on 
the quantitative information derived from it. 
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Figure 3 


Block diagram of optical processing and pattern recognition system 





PATTERN 
ORIENTATION 
& SPATIAL 
FREQUENCY 
INFORMATION 


SPATIAL FREQ. (RING 57 


ORIENTATION (WEDGE) DATA 



Figure 4 . 


Detail of portion of Fig. 3 enclosed 
in dashed lines. 
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3.2 


Procedures and Measurements 


3.2.1 Optical Processing System Description 

In this section we will present a description of the optical processing system 
used for spatial frequency analysis of the ERTS images. The optical processor has 
three main elements: a laser, optics, and a Recognition Systems, Inc,, Diffraction 
Pattern Sampling Unit (DPSU). The system configuration is shown in Figure 5. An 
ERTS— 1 70 mm positive transparency is used as the input for this system. The optical 
processing system can be regarded as a two step system. First, an area of the ERTS 
transparency (sample area) is illuminated by the incident laser beam. This beam is 
focused by the lens so that the point source produced by the spatial filter is imaged 
at a distance z + f in front of the lens. The light distribution in this plane is pro" 
portional to the Fourier transform of the transmittance function which describes the 
portion of the ERTS image illuminated by the beam. Second, the intensity distribution 
(spatial spectrum) of the ERTS image is sampled by the DPSU. 

The DPSU consists of a 64 element photodiode array (shown in Figure 6) 
used to detect the light intensity incident upon each element, and electronics which 
amplify and digitize the output from each diode in the array. The diode array is 
composed of 32 wedge-shaped photodiodes and 32 annualar ring photodiodes. The 
average intensity for each photodiode is recorded. These data are then used in a 
computer program and are calibrated, printed, and plotted. 

The spatial frequency in the transform plane is related to other system 
parameters by: 

s = r/dA (2) 

where s = spatial frequency at a point r in the transform plane 
r = distance in transform plane measured from the optical 
axis 

d = distance from image transparency to detector 
^ - wavelength of laser radiation = 6328 A 
The spatial frequency from this calculation is converted to ground spatial frequency 
using image to ground scale. The resulting curves which are plotted by the computer 
program are then intensity vs. frequency and intensity vs. angle. These are plotted 
in terms of ground spatial frequency in cycles per kilometer and direction in compass 
degrees from north. 
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In order to calibrate the detector array, the array was illuminated with a 
uniform intensity, plane wave source,. The intensity of the illumination was recorded 
for each element in the array and the measurements were used as the basis for the as - 
signment of an area/sensitivity calibration factor to each element. 

The mathematical description of this system is given in Appendix A. 
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3*2,2 Data Taking Procedure 

The procedure used to record the data for each sample area Is outlined be- 
low; 

(1) Set DPSU digital readout to zero, (Adjust 2nd Stage Balance Control ~ 
ring 30 switch on- set display = .000 at 1000 X gain.) 

(2) Place image on image holder, (image is aligned with respect to 
holder) 

(3) Choose sample area by moving ERTS image in image (x,y) plane. 

(The area on the ERTS image is chosen initially with respect to an 
x,y coordinate system.) 

(4) Align detector with respect to transform. (This is accomplished by 
moving the detector array in the transform plane with the center 
element (ring 1) switched into the circuit until the maximum reading 
is obtained on the digital display. This is a very critical adjustment 
that aligns the detector array with the optical axis of the system.) 

(5) Record data by switching each element of the array, in turn, into 
the circuit. 

(6) Enter data onto computer cards. 

(7) Run data calibration program. (Output Is printed and automatically 
plotted.) 
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3.3 


Experimental Results 


The spatial frequency curves obtained for each sample area were modified 
to enhance their ability to display features of interest. An average spatial frequency 
curve was calculated for the 80 non-snow sample areas. The point-by-point ratio 
between this curve and each non~snow sample area frequency curve was determined 
and these values were plotted. A similar procedure was used for the snow sample 
area frequency curves. The relative amplitudes of the various frequency components 
are greatly emphasized using this method and it will be shown in the next chapter 
that these modified spatial frequency curves characterize the sample areas very well. 

At this point in the analysis of the sample areas each sample has been des- 
cribed by a modified spatial frequency curve and an orientational curve* The steps 
involved in going from the ERTS sample areas to the modified spatial frequency and 
orientational curves are illustrated in the following section. 
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3.3.1 Analysis of Results for Sample Area Example 

This section is intended to demonstrate how specific features In a sample 
area can be related to the spatial frequency and orientational information derived for 
it. Figure 7(a) shows a portion of ERTS-1 image no. T076 - I6393~5 , The circular 
area is sample area F“7 from the Flint Hills region of Kansas. 7(b) shows the dif“ 
fraction pattern associated with the sample area shown in (a) as obtained with the 
optical processor. This pattern can be related to features in the sample area. The 
main feature of the diffraction pattern is the very bright distribution (labelled as A) 
extending vertically through the pattern. This distribution is due to the scan lines 
in the ERTS MSS images associated with the manner in which these images were 
obtained. Distributions which extend vertically in this manner are not detected 
with the DPSU, however, since a 20° dead space (see Fig. 6) exists in this region. 

It should be noted that because of the nature of the diffraction pattern, scan lines, 
which appear horizontally across the image, produce this vertical distribution. The 
intensity distribution in (b) is sampled using the DPSU described in the last section. 

The normalized, calibrated spatial frequency (1) and orientational (2) curves 
are shown In Figure 7(c). As noted, these plots are obtained by recording the light 
Intensity at each photodiode in the detector array and then using these data In a 
computer program. The spatial frequency curve, obtained from the ri ng~shaped 
photodiodes, gives the intensity of various spatial frequency components corresponding 
to features In the image sample area. The normalized intensity (shown here on a 
logarithmic scale) is recorded versus ground spatial frequency in cycles/kilometer. 

The orientational curve, obtained from the wedgeshaped photodiodes, indicate the 
orientation of features in the sample area. The normalized (with respect to the 
largest value) intensity is recorded versus the orientation clockwise in degrees from 
image north. 

Figure 7(d) shows the modified spatial frequency curve obtained from the 
spatial frequency curve in (c). The amplitudes of the various frequency components 
are greatly emphasized in this modified frequency curve. The curve in (d) shows that 
frequency component values around .4 and 1 .0 cycles/km tend to predominate. There 
are also peaks in the spatial frequency curve shown in (c) around .4 and 1.0 cycles/km. 
We would suspect that these peaks are caused by characteristic features in the sample 
area with these frequencies. The sample area in (a) shows that this is true. A band 
which is light in tone is apparent in this image. This band represents the Kansas River 
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Figure 7. Example Showing Progression of Steps From Portion of ERTS Image No. 

1076-16393-5 (a); to Diffraction Pattern (b); to Digitized Computer Plots 
(c); to Modified Spatial Frequency Curve (d). 
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flood plain which is an area of intensive agriculture. This band is approximately 
2.5 kilometers wide and hence it has characteristic frequencies around 1/2.5 = 0.4 
cycles/km. The distribution around 1.0 is apparently produced by characteristic 
stream pattern frequencies. The orientational curve, curve (2) in 7 (c), also has 
a number of features which can be related to features in the sample area. 

The two dominant peaks in this curve are labeled B and C. The peak at 
B, at approximately 35°, appears to result from the orientation of the Kansas River 
flood plain. This plain is oriented at approximately 30° to 35° clockwise from west 
in the sample area. Thus, because of the nature of the system, this plain produces 
the distribution at B in (b) which is oriented at 90° away. The peak at C in curve 
2, 7(c) which occurs centered around 90° is due to the northward orientation of the 
field patterns and the general orientation of the major streams in the sample area. 

The distribution at C in (c) is also relatively wide indicating the variation in stream 
orientations around 90°. This example should illustrate how features in each sample 
area can be related to features in the diffraction pattern and in turn related to features 
in the spatial frequency and orientational curves. 
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3,3,2 Example of Results for the Eight Categories 

Several of the ERTS sample areas are shown in Figures 8(a) through 13 (a). 

The resultant spatial frequency and orientationalcurves are shown In (b) for these areas, 
and the modified spatial frequency curves are shown in (c). (Peaks in the spatial fre“ 
quency curves in (b) around ,6 to .8 cycles/kilometer are usually caused by multiple 
scatter in ring 13 of the DPSUand hence are unrelated to image features,) 

Figures 8(a) and 9(a) show two sample areas from the Flint Hills region of 
Kansas, Stream patterns In these two sample areas are well expressed and result in 
high peaks in the orientational curves in 8(b) and 9(b) corresponding to the orienta" 
tion of the major streams. The modified frequency curves in 8(c) and 9(c) demonstrate 
characteristics of the modified frequency curves from the Flint Hills, the Blue Hills, 
and the Smoky Hills, The slope of these curves in the region beyond approximately 
1,0 cycles/kilometer is positive and demonstrates the relatively larger high frequency 
content for sample areas from these categories. 

Figures 10(a) and 11(a) show sample areas from the Osage Plains region 
where Figure 11(a) is a sample area with snow cover. There are numerous peaks in 
the orientational curves in 10(b) and 11(b) which is characteristic of orientational 
curves from sample areas from the Osage Plains, The modified frequency curves for 
both sample areas 10(a) and 11(a) fall off at higher frequencies in contrast to the 
curves in 8 (c) and 9 (c). 

Figure 12 (a) shows sample area G~5 from the northeastern part of Kansas. 

The orientational curve in 12(b) has high peaks at 90° and 180° indicating the dominant 
north-south, east-west orientation of the section roads and agricultural patterns in 
this area. Figure 13(a) shows the same sample area as in 12(a) with snow cover. The 
snow cover emphasizes the stream patterns in this case and the orientational curve, 

13(b), has high peaks around 50° -75° corresponding to the major stream orientations. 

The modified frequency curves shown in 12(c) and 13(c) again fall off at high fre- 
quencies. 

These examples are intended to show how features in the spatial frequency 
and orientational curves can be used to categorize the sample areas. In chapter 4 
it will be shown that features such as those shown here can be used to discriminate 
between various physiographic categories. 
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Figure 8. (a)Sample Area F-4 ; (b) Optical Processor Curves ; (c) Modified Frequency Curve 
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Figure 9. (a)Sample Area F-6 ; (b) Optical Processor Curves; (c) Modified Frequency Curve. 
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Figure 11. (a) Sample Area 0—20; (b) Optical Processor Curves; (c) Modified. Frequency Curve. 
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Figure 12. (a) Sample Area G-5; (b)Optical Processor Curves; (c)Modified Frequency Curve. 
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Figure 13. (a) Sample Area G—15; (b)Optical Processor Curves; (c)Modified Frequency Curve. 



3.3,3 Multispectral Example 

Figure 14 shows portions of ERTS-1 images 1006~16502“4,5,6,7 and the 
sample areas M“4, M _ 5, M“6 and M - 7. The spatial frequency and orientational 
curves obtained with the optical data processing system are shown next to the re- 
spective samples. These images are from MSS bands 4(.5-.6 micrometers; 5(»6 - .7 
micrometers); 6(.7-.8 micrometers); and 7(. 8-1.1 micrometers). The sample area 
shown here is in northwestern Kansas and contains a parallel drainage system which 
trends NW-SE. These four bands were analyzed in order to obtain quantitative in- 
formation on the way in which features, specifically drainage patterns, are expressed 
in the various bands. 

The expression of the drainage pattern in the four images for the four bands 
is a function of the Iand - use in this sample area from the High Plains region of Kansas. 

The rough stream valleys here are often used as pasture land and the level uplands are 
used for cultivation. The main feature that concerns us here is the difference in 
expression in the images between the grassy pastureland and the cultivated ground 
in the uplands. An idea of the relationship between the expression of the land-use 
categories in the images and the actual land-use can be obtained from curves of 
relative spectral radiance for various agricultural scenes. Curves of this type 
are shown in Remote Sensing (1970). Specifically, Figure 7,page 310shows howthe relative 
spectral radiance of four types of agricultural scenes changes as a function of wave- 
length. The main difference here is between the radiance of the various scenes in 
the visible and the infrared part of the spectrum. The curves shown in this figure 
demonstrate that the radiance of the grassy pastureland should be greater in the in- 
frared bands than in the visible band and conversely, the radiance of the bare ground 
should be greater in the visible than in the infrared. This relationship is apparent 
In the image specifically where the drainage patterns are dark in tone in the visible 
bands (4 and 5), and are displayed as light in tone in the two infrared bands (6 and 
7). Image features as they are displayed in the spatial frequency and orientational 
curves are described in the next paragraph. 

The orientational curve for the band 4 image contains a peak at 100° which 
is due to the field patterns apparent in the image. The drainage pattern in this image 
is not enhanced as much as in the band 5 image due to the weaker reflectance of the 
field patterns in band 4. Because of the number of field patterns in the sample area. 
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Figure 14 Portions of ERTS— 1 Images 1006“ 16502-4,5,6, 7 and Optical 
Diffraction Curves for Sample M without Snow. 
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the contribution due to the field patterns slightly dominates in the orientational curve. 
The parallel drainage pattern produces the secondary peaks at approximately 75° and 
85°. The drainage pattern in the band 5 image is very apparent due to the enhance - 
ment effect provided by the field patterns. Correspondingly, the orientational curve 
displays the peak at 75° due to the drainage direction , and displays a secondary 
peak associated with the field patterns at 100°. The band 6 and band 7 images are 
similar and produce similar results. Both the drainage patterns and the field patterns 
appear to be subdued in these two images apparently because the medium reflectance 
of the fields again does not enhance the drainage patterns. Again because of the 
number of fields in the sample area, however, the orientational curves both have 
sharp peaks at 100° that result from these field patterns. The drainage pattern for 
both of these images does not contribute significantly to the orientational informa- 
tion displayed in the curve. The spatial frequency curves for all four images are simi- 
lar. However, since features apparent in the sample area are dilineated more sharply 
in the band 4 and band 5 images, the spatial frequency curves for these two bands 
do display a slightly greater high frequency content. 

This example gives a quantitative view of why, in order to emphasize geologic 
ground features, band 5 images were chosen for analysis in this investigation. 
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4.0 CLASSIFICATION RESULTS AND ANALYSIS 

4.1 Derivation of Pattern Descriptors 

To reduce the amount of information present in the modified spatial fre- 
quency and orientational curves, various data processing schemes were developed 
to extract parameters from them. These parameters describe the geologic or physio - 
graphic features giving rise to fluctuations in the curves. To determine how well 
these parameters characterize each sample area, they may be used to categorize 
the sample areas by using the parameters as input to pattern classification algorithms. 

The range of frequencies was divided into two bands. Band 1 contains 
spatial frequencies between 0.0 and 0.9 cycles/km and band 2 contains spatial fre- 
quencies between 1.1 and 2.8 cycles/km. This division essentially separates in- 
formation due to high frequency fluctuations such as occur in stream patterns in rough 
terrain, and low frequency information such as occur due to field patterns. This 
division of frequencies appears to correspond to a natural break at approximately 
1.0 cycles/km. 

Similarly, the range of orientational data was divided into four sectors — 
each corresponding to 40 degrees (see Figure 6). Sector 1 provides data on pattern 
orientations that produce distributions between 25 and 65 degrees clockwise from 
north; sector 2: 65 - I05 degrees; sector 3: 105-145 degrees; and sector 4: 145-185 
degrees. 

Parameters are extracted from the spatial frequency and orientational curves 
in each range . 

The following parameters were calculated from the spatial frequency and 
orientational curves in each of the bands or sectors listed above. 

Orientational Curve Parameters: 

DAV Average value of the curve in the 40° sector 

RMSD Root Mean Square Deviation of curve from overall average in each 

sector 

AREA A. Area of curve above overall average in sector 

B, Area of curve below overall average in sector 
PEAK Number of peaks or 'spikes' in curve in sector 

PAKS Number of data points above average in sector 
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Spatial Frequency Curve Parameters: 

(These parameters are derived from the modified spatial frequency curve) 


DAV 

Average value of the curve in each frequency band 


RMSD 

Root Mean Square Deviation of curve from overall average 

in each band 

AREA 

A. Area of curve above overall average in band 

B. Area of curve below overall average in band 


DARA 

A. Area of curve above the value 1.0 in band 

B. Area of curve below the value 1.0 in band 


DYNR 

DYNamic Range of curve in band 


SLOPE 2 

The slope of the regression line calculated for the curve in 

band 2 


Each parameter for all 80 non-snow and 68 snow sample areas was plotted 
versus its appropriate sample area category as obtained from Figure 1 . Three of these 
parameter plots are shown in Figures 15, 16, and 17, Two separate parameters are 
listed and one of these is repeated for the 68 snow sample areas. The notation 'R 1 
or 'W' preceding the parameter name denotes the parameter as a spatial frequency 
(ring) or orientational (wedge) parameter respectively. 
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PHYSIOGRAPHIC CATEGORY 


Figure 15. Parameter Plot Using the Parameter (R) SLOPE2 




(R) SLOPE2 - AMPLITUDE (SNOW) 



PHYSIOGRAPHIC CATEGORY 


Figure 16. Parameter Plot Using the Parameter (RJSLOPE2 for the Snow-covered Sample Areas. 
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4.2 Physiographic Patterns Identified 

Decision boundaries may be inserted into the parameter plots shown in the 
last section based on the distribution of the amplitudes of the parameter. If there is 
a clear relationship between the amplitude of the parameter and the physiographic 
characteristics of the sample area, then this decision boundary will discriminate be - 
tween 2 large-scale physiographic ground patterns. 

This was done for the parameter labelled SLOPE 2 and the results of this 
scheme are shown in Figure 18. As mentioned previously, SLOPE 2 is merely the 
slope of the regression line for the modified spatial frequency curve for the frequen- 
cies in band 2. As can be seen from the plot showing the parameter SLOPE 2 (Figure 15), 
a decision boundary may be inserted at approximately zero slope (+.001) which 
effectively separates the various categories. Sample areas from categories of re- 
latively large high frequency content (Flint Hills, Red Hills, Smoky Hills) yield a 
positive slope, whereas sample areas from the other categories yield a negative 
value of this parameter. 

This identification experiment resulted in an accuracy of 92.5% with only 
6 incorrect identifications based on the original categorization of the sample areas. 

Sample sites represented by a square were identified as belonging to the category 
comprising the Flint Hills, Red Hills, and Smoky Hills; sample sites represented by 
a triangle were identified as belonging to the category comprising the Arkansas 
River Valley Lowlands, High Plains, Osage Plains, Blue Hills, and the Glaciated 
Region. 

Figure 19 illustrates how a further classification based upon the original 
physiographic regions can be made. The figure shows a scattergram made by plotting 
values of the parameters DARAN 2 versus values of SLOPE 2. Values on this graph 
corresponding to sample areas in the Red Hills, Flint Hills, and Smoky Hills are 
designated by an open circle; values corresponding to sample areas in the Blue Hills 
are designated by a black circle; values corresponding to sample areas in the Glaciated 
Region and Osage Plains are designated by an open square; and values corresponding 
to sample areas in the High Plains and the Arkansas River Valley Lowlands Regions 
of Kansas are designated by a black square. 

The categories listed above each appear to be distributed in separate por- 
tions of the graph. It is apparent that we may insert decision boundaries into this 
graph to separate the various categories as was done in the one -dimensional case 
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kilometers 

Figure 18. Classification of Sample Areas in Terms of a Decision Made on the Value of the Parameter SLOPE2 
for Each Sample Area. 





shown previously. If we now define a region for each category in terms of these 
boundaries, we may determine how accurately these two parameters classify each 
sample area. This was done and the results of this classification are shown in Table 1. 
The overall accuracy in terms of a correct identification for each sample area is 
80% for this case . 

The sample areas are shown in Figure 20 with respect to their identification 
based on Figure 19. Figure 20 shows that although the sample areas were identified 
with an accuracy of only 80%, the sample areas appear to be clustered in Kansas 
in terms of their identified categories. For instance, A-4, A -5, and A“6 are classi- 
fied with the Blue Hills to the north, and H“8 is classified with the Blue Hills im- 
mediately to the east of it. The category designation in Figure 20 is the same as 
in Figure 19. 

Various algorithms exist to effect the classification made here automatically. 
One such scheme given in Fukunaga (1972) uses a regression type algorithm to assign 
samples to various categories. A program based on this algorithm and developed by 
Haralick and Shanmugan (1973) was used in an attempt to classify the sample areas, 
but we had only limited success with this program because of the relatively small 
number of samples used here. 
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TABLE 1 


Result of Classification of Sample Areas Based on Figure 19, 


CATEGORY 

HIT 

MISSED 

TOTAL 

SAMPLES 

% 

ACCURACY 

R,S,F 

25 

5 

30 

83 

B 

7 

3 

10 

70 

G , 0 

19 

1 

20 

95 

H, A 

13 

7 

20 

65 

TOTAL 

64 

16 

80 

80 
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Figure 20. Classification of Sample Areas Based on Figure 19. 






4.3 


Information Content of Frequency Bands 


In an effort to determine the bond of spatial frequencies that contain the 
most information with respect to discriminating between physiographic categories, 
an additional experiment was performed. Data from the optical processing of some 
of the ERTS sample areas were obtained using a reduced spatial frequency scale. 

The scale in the transform plane was changed in order to obtain spatial frequency 
information on frequencies between 2.8 and 5.8 cycles/km. 

Two of the reduced scale modified spatial frequency curves obtained in 
this experiment are shown in Figure 21 (a), and (b). These curves are from the same 
sample areas shown in Figures 9, and 12. The three frequency bands over which 
parameters are calculated are also illustrated in this figure. Band 1 contains fre- 
quencies between 0,0 and 0.9 cycles/km; band 2 contains frequencies between 
1 .1 and 2.8 cycles/km; and band 3 contains frequencies between 3.1 and 5.8 cycles/ 
km. The feature parameters described previously were obtained for the additional 
frequency band (band 3). The parameters calculated for band 3 appeared to provide 
approximately the same information with respect to discriminating between the physio- 
graphic categories for some of the parameters. In general band 3 did not work as 
well as band 2 in separating the categories. 


38 



RELATIVE INTENSITY 



Figure 21. Examples of the Modified Spatial Frequency Curves Showing the Curves 
in All Three Frequency Bands for Sample Areas G-5, F-6. 





5.0 


CONCLUSIONS 


The evidence presented in this paper establishes that optical data processing 
of ERTS-1 images can be used very successfully in identifying large-scale ground pat- 
terns in Kansas. An Important point is that the sample areas that were used in this 
investigation are not uniform. That Is they do not display merely one dominant 
characteristic. Because each sample area is approximately 37 kilometers in diameter, 
the terrain may vary substantially within each sample area. Even given this 'noise 1 
associated with the sample areas, this method works very well to Identify the sample 
areas in terms of the original physiographic categories. 

The manual interpretation of the spatial frequency and orientational curves 
shown here for one sample area and presented for all the sample areas in Ulaby et al. 
(1973), demonstrated that this method of analysis can provide an unbiased means of 
evaluating the orientations and natures of stream patterns and other features apparent 
in the Image . 

The band of frequencies between 1.1 and 2,8 cycles/km appears to contain 
most of the information useful in discriminating between the various physiographic 
categories. The nature of the information contained in frequencies below 1.1 cycles/ 
km is not apparent with respect to classifying the categories; the nature of the in- 
formation in the frequencies from 3.1 to 5.8 cycles/km appears to discriminate be- 
tween categories almost as well as the band from 1.1 to 2.8 cycles/km but only for 
certain parameters. Sample pattern orientational information is not as useful in dis- 
criminating between categories as the spatial frequency information. 

Snow cover suppresses patterns due to agriculture and vegetation while 
enhancing the topography. Such topograhic enhancement aids manual geologic 
interpretation and stream pattern analysis. However, the processing of the sample 
areas with snow-cover did not provide any help in terms of identifying the sample 
areas. In fact, snow cover seemed to obscure the information useful in discrimina- 
ting between categories. It is possible that the snow cover obscured the contribu- 
tions due to the high frequency information in the band from 1 .2 to 2.8 cyclesAm. 

It is important to note here that the method of analysis used in this investi- 
gation is very general* The same methods and algorithms could be applied in general 
to various other image analysis or categorization problems. The method of analysis 
shown here is image independent and could be completely automated. 
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APPENDIX A. Optical Processor Equations 


The Fourier transforming property of a simple lens 
system was demonstrated in section 2.1. Now we consider 
the more general case used in this investigation and il- 
lustrated in Figure 5 and again here in Figure A-l, where 
a point source at 0 is imaged by a single positive lens at 
a point O’ in a plane at U^. 


U ,4 U ? U 2 U .l 

d »-] |*D-*| 



Figure A-l. Configuration Assumed for the Derivation of 
Equation A-6. 

The transparency in this case is placed behind the lens. 

We will show here that this general configuration also 
demonstrates the Fourier transforming property and illustrate 
the origin of Equation 1. 

The complex distribution in front of the lens in 
plane using a paraxial approximation to a spherical wave, 

is : 

« 5 y t ) = A e*p(jkd 0 ) expl + y, 1 )) (A-l) 

The distribution behind the lens in plane given the 
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distribution in is (from eqs. 5-2, 5-10, pp. 78-81 in 
Goodman) : 

^2^2)72.) — A exp (j k<J 0 ) e *p£jfJ~ ( * t Z + y t Z )J 

X <rxp(jkr\t>) exp L- (? 2 2 + 72 *')] (a- 2 ) 

= A e%p CjkA B ) e-^pCjknb)eypl-j-j-(~f~^- )(* z ^ 7-2 )] 

where x 2 = x j ; y 2 “ y^; from the thin lens approximation. 

Now let: Aexp(jkd ) exp(jknD) = A 1 ; and use the lens law 

J L = -L = . J _ 

■f. . Jo 'fi f, +*, 

to obtain: 

V x (** >72) = ^ 'e^pr-J^C ^ 2 + (A_3) 

This equation represents the complex distribution behind 
a lens of focal length illuminated by a normally in- 
cident plane wave of amplitude A 1 -. 

At plane / the transparency, represented by a com- 
plex transmittance function t , is multiplied by the in- 
cident distribution to obtain (Goodman p. 88) : 

U 3 b 3 ^s) = U^y 3 )A^ 

X erp [~j£j ( x 3* + y3 Z )J < A_4 ) 

where the distribution is scaled appropriately and the 
f 2 f 2 

factor P(— x^, — ^p-y^} represents the area on the image 
transparency illuminated due to the projection of the pupil 
function of the lens down the cone of light rays onto the 
transparency . 

If Fresnel diffraction is assumed from the object 
plane to U^, then (Goodman p. 60): 
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\j, o„ , * ¥ ) = V 4 V ) 

^ J/ ^(*3 )^ ex P )] 

- 00 

X ^ V Y ^3 3 ^*.? (A_5) 

Finally, substituting eq. (A-4) into eq. (A-5) : 


U*(V*,7y) = 2— exp[ j^O/ if 

*o 

^ ff t a (x 3 ,y 3 ) P(-^X 3 , 

-OO 


^ expT-jf^CXjXv +-y 3 ^y)]c /^ 3 c/-^ (A _ 6) 

where a constant phase factor has been dropped. 

If the quadratic phase factor is neglected, eq. (A-6) 
states that the distribution in plane is merely the 
Fourier transform of the trasmission function describing 
that portion of the transparency illuminated by the pro- 
jected light rays and scaled appropriately. Thus, this 
general lens system also demonstrates an inherent Fourier 
transforming property. Equation 2 is obtained from eq. 

(A-6) by noting that this Fourier transform has a scale 


factor associated with it based on the geometry of the 

system. That is, the transform must be evaluated at fre- 

X 4 y 4 

quencies f - — ~r , f = — — -=■ , . . _ , 

x Xd y Xd , or combining these we obtain 

r 

s = which is equation 2, where s equals the spatial 

frequency at a radius r from the optical axis in the 


transform plane. 
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APPENDIX B. 


Equipment 



Figure 33-1. Optical Processing Equipment. 


1 . 

Laser 


Spectra-Physics Model 134 HeNe 
laser. Output power--4mw. 
Wavelength — 6328 Angstroms. 

2. 

Microscope 

objective 

Acts as spatial filter and beam 
-expander. Output is an Airy pat- 

3. 

25 micron 

pinhole _ 

tern. 

4. 

Lens 


The lens (focal length = 16.8 cm) 
images the point source formed at 


3 in the plane at 6. The lens 
iris diaphragm is adjusted to pass 
the central lobe of the Airy pat- 
tern formed at 3. 

5. Image holder The image is inserted and moved 

in this plane to select the de- 
sired sample area. 

6,7. DPSU Recognition Systems Incorporated 

Diffraction Pattern Sampling Unit 
(DPSU). 64 element photodiode 
array at 6, (see Figure 6), samples 
the diffraction pattern formed in 
this plane. The unit at 7 ampli- 
fies the output from each photo- 
diode. Each photodiode can be 
switched into the circuit and the 
resultant amplitude displayed. 
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TABLE B— 1 


Average Radii of Photodetector Annular Rings (see Fig. 6) 


ring no. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 


r (mils) 


0.0 

8.0 

14.3 

21 . 8 

30.3 
40.0 
50.6 
62.2 
74.8 

88.3 
103 . 
118 . 
135 . 
152 . 
170 . 
189 . 
209 . 
229 . 
251 . 
273 . 
296 . 
321 . 
346 . 
372 . 
398 . 
426 . 
454 . 
484 . 
514 . 
545 . 
577 . 
610 . 


r - Average radius from center 
of detector. 


(These are the average dis- 
tances from the center of the 
detector array for each ele- 
ment in thousandths of an inch 
as calculated from the manu- 
facturers data.) 
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( From Ulaby et al. (1973) ) 


Kansas lies within the Stable Inferior, a large geologic province of North 
America occupying most of the region between the Appalachians end the Rocky 
Mountains. This area has suffered little in the way of intense tectonic activity 
since early Cambrian sediments were first deposited approximately 600 million years 
ogo. The area which is now Kansas was the site of shallow seas during a good 
portion of this period. The result has been the formation and preservation of 
sedimentary formations (sandstones, shales, and limestones) which cover the much 
older pre-Cambrian basement rocks. The total thickness of these sedimentary rocks 
is nowhere more than 9500* . This sedimentary cover Is thin when compared to 
sedimentary thicknesses in other parts of the country where sedimentary basins may 
be 20,000* to 30,000* deep. 

The sedimentary formations in Kansas are relatively thin and exhibit a slight 
but persistent westward dip at the surface. The present day streams and rivers of the 
state flow in a generally eastward direction from the piedmont of the Rocky Mountains 
to the Missouri and Mississippi Rivers. This reflects the uplift and eastward tilting 
of much of the state which occurred with the formation of the Rocky Mountains. 

Thus the surface elevation of the state increases from the eastern border to the 
western boundary ranging from less than 800' to more than 4,000* . In traveling 
In a westward direction across the state, one not only gains elevation but crosses 
the outcrops of progressively younger sedimentary rocks. Because of the westward 
dip of the rock units and the subsequent erosion, they are arranged in a stairstep 
fashion with each upward step representing a younger sedimentary unit. 

The lithologic makeup of a rock unit, its structural attitude and the 
weathering forces operating upon it determine the landforms that will form in the 
area of its outcrop. Because of the westward dip of sedimentary rocks in Kansas, 
and their layer-cake arrangement, much of the surface of the state is characterized 
by a series of eastward facing escarpments or cuestas. These features occur where 
a relatively resistant unit (such as a limestone) crops out. Such a unit tends to 
produce a londform having a steep eastern face end gentle western slope. This 
gentle back slope often extends across the outcrop of an overlying, less resistant unit 
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(such as a shale) to the base of the escarpment formed by the next resistant unit. 

These esccrpments are dissected by the erosive effects of the many eastward flowing 
streams and as a result have an irregular appearance on a map or aerial photo. But 
On a larger scale they maintain a general north-south orientation which is at right 
angles to the prevailing dip direction of the exposed rock units. 

Lithologic differences are responsible for the various physiographic regions 
generally recognized in Kansas which are shown in Figure 1 . In the eastern part of 
the state are the Osage Plains which are made up of a series of ecstward~facing 
escarpments formed by outcropsof Pennsylvanian limestones and shales. The Flint 
Hills are actually a large escarpment formed by outcrops of a series of chert-bearing 
Permian limestones that are very resistant to erosion. The Dakota sandstone outcrops 
(Cretaceous) form the Smoky Hills upland in the North Centra! part of the state . 
Likewise, the Blue Hills are formed by outcrops of Upper Cretaceous limestones and 
chalks. 

Beyond the Blue Hills lie the High Plains which are formed by the accumulation 
of sediments derived from the erosion of the Rocky Mountains to the west. Numerous 
aggrading streams swept eastward during the Tertiary carrying and depositing sand 
and gravel and forming a vest oufv/ash plain, that is toacy the present land surfece. 
Even younger deposits occur in the various prairies and lowlands associated with the 
Arkansas River. Some of these areas are covered with wind blown sand in the form 
of dunes, both stabilized and active. In the south-central part of the state are the 
Red Hills or Cimarron breaks marking the border of the High Plains in the vicinity 
of the Cimarron River, which together with Its tributaries eroded Into the red Permicn 
siltstones and shales which underlie the High Plains in that area. The extreme 
northeast corner of the state was occupied by the Kansan glccier during the 
Pliestocene. As a result, the landscape was resculptured to some extent and the area 
was covered to varying depths by glacial deposits. 


STRUCTURAL GEOLOGY 


The structure of Kansas Is subtle for the most part, seldom being dramatically 
expressed on the surface. This should be kept in mind when considering remote 
sensors as tools in mopping the geology or structure of Kansas. 
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Although faulting and folding are not intense or widespread in Kansas, joints 
are. According to Billings (1942), “Joints may be defined as divisional planes 
or surfaces that divide rocks and along which there has been no visible movement 
parallel to the plane or surface". As Merriam (1963, p. 254) states, * Little work has 
been done on jointing in Kansas, although the joints are extensively developed . 

Word (K.G .5. Bulletin 191, pt. 2) in 1968 did a study of joint patterns in the 
Southern Flint Hills and he concluded, among other things, that the joint patterns 
measured showed a close correlation to regional tilting and may have been produced 
at the time of the tilting. He also states that, "The present drainage patterns appear 
to be closely related to and may be determined by the joint pattern” , (Ward 1968, 
p. 21). He adds, like Merriam before him that "more work concerning midcontinent 
joint systems is justified" (Ward 1968, p. 21). 

Stream patterns have long been known by geologists to reflect the underlying 
structures, faulting, folding end jointing. Such deformation tends to rupture 
competent formations creating planes of weakness along which weathering activities 
are occelerated and through which streams have a tendency to flow, taking advantage 
of the destructive work already cone for them. Thus, by studying topographic parterns 

in an area, insight can be gained concerning the geologic structure. 

Many studies have been performed by geologist using aerial photographs 
as an aid in structural analysis. Kelley (1960) mapped regional fracture systems 
for a large area of the Colorado Plateau using aerial photography. Boyer and 
McOwen (1964) working in Texas established a relationship between fracture 
patterns observed on the ground and linear features on aerial photographs. Likewise, 
in the Appalachian Plateau, Lattman (1958) established a correlation between 
bed rock joint systems and linear features on aerial photographs. These studies 
and many others like them involve the visual detection, measurements, and 
evaluation of linear features or lineaments. Such studies are limited by the 
interpreter's ability to detect lineaments and may be handicapped by a biased 
evaluation of their significant. Spatial frequency analysis may provide a means 
of detecting and measuring linear features on imagery which would not ordinary be 
detected by visual means. In addition, such analysis would not be guided by any 
preordained knowledge of the geologic structure, and the results would be unbiased. 
The small scale of the imagery used in this study negates the detection of actual 
fractures or joints on the ground. However, the layer linear features associated 
with stream traces and topographic alignments can be detected and will provide 
most of t ho information concerning geologic structure. 
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LAND USE 


Lond”Use In Kansas Is predominately devoted te agriculture. However / the 
type of agriculture practiced varies across the state, due to climate, soil, landform, 
and availability of water. Agricultural land - use can be correlated fairly well with 
the physiographic regions of Figure 1. The area east of the Flint Hills, namely the 
Osage Plains and the Glaciated Region, is characterized by mixed farming of the 
corn belt type, with generally small fields and pastures and a variety of crops grown 
including corn, soybeans, milo, etc. Hay is an important crop in some areas of the 
Southern part of this area and stock farming is a common practice . To the west lies 
the Flint Hills region with Its areas of bluestem prairie . It is predominately a 
ranching area specializing In finishing out transient cattle brought in from western 

end southwestern ranges. The Blue Hills and Smoky Hills are also large cattle grazing 
regions, especially in the rougher areas where cultivation is impractical. In this 
some regard the Red Hills region of the south-central part of the state is also an 
important rcnching area. Much of the arable area in the western two-thirds of the 
state is devoted to the raising of wheat, with the most extensive whect growing 
areas in the level and fertile lowlands associated with the Arkansas River. Much 
wheat is also grown on the level uplands of the High Plains and in favorable areas 
in the Smoky Hills and Blue Hills. However, dry farming is in common practice 
in the western part of Kansas, where rainfall is deficient. Dry farming involves 
practices designed to catch and conserve the available moisture. Toward this end, 
dry farming often involves the fallowing of fields for one or two yeers in order 
to build up a reserve of soil moisture. Thus in any given year in the western pert 
of the state, a sizcble portion of the cultivated land will be free of planted crops. 

Much of the High Plains is also devoted to the grazing of cattle especially the rougher 
lands along streems. The southern High Plains of Kansas is an important grain 
sorghum growing area. 

Irrigation is important in ports of western Kansas . The largest and most 
extensive area of irrigation is centered around Garden City, reaching from Scott 
City southward into Meade County and westward along the Arkansas River. This 
area is underlain by a sedimentary basin containing thick deposits of the Tertiary 
Ogallala formation. This formation is largely sand and gravels derived from the 
Rocky Mountain in Tertiary time. Its porous end permeable nature moke it an excellent 
aquifer (water yielding formation) and it feeds the many irrigation wells in this 
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orea. Several crops are grown here including com, wheat, soybeans, alfalfa, and 
sugor beets. Other irrigated areas in the High Plains regions also rely on the ground 
water stored in the Ogaflala, The more recent deposits of alluvial materia! in 
the valley of the Arkansas River are also important aquifers. The tell-tale signs 
of irrigation on areal photos are the circular fields produced by pivitol sprinkler 
systems. The circular fields are quite large (1/2 mile in diameter) and are discernibl 
on ERT5 imagery as well. 

The Arkansas River valley also contains areas of sand hills. In most of these 
areas, the sand hills or dunes are stabilized and covered by natural vegetation. 

Some areas are active however with dunes in formation and smaller areas of wind 
erosion colled blowouts. Because of the rolling topography of these sand hill 
tracts, many are uncultivated and used as grazing areas. 

On an area basis, non-agricuitural land use is, of course, secondary in 
Kansas. The remaining land area is largely tied up in cities and towns, reservoirs 
end their surrounding management areas, military installations, wildlife refuges 
and mining areas. These land uses are largely self-explanatory . However, land 
use related to mining activity requires further elaboration. 

Among the minerals and rocks mined In the state are lecd, zinc, coal 
gypsum, salt, volcanic ash, limestone, sand, gravel, and clay. In addition, the 
state is an important oil and gas producer. Of these activities; salt, lead, and 
zinc mining is performed under ground. Most of the other products are mined by 
quarrying operations that are generally small In areal extent. The exception is 
the procurement of coal which is done by strip mining. Large areas in southeast 
Kansas bear the effects of strip mining in the form of long parallel mounds of dirt 
and rock which represent the overburden removed to reach the underlying coal 
seam. Lakes occupy many of these abandoned strip pits today and are used as a 
recreational resource in the area. The reclamation of strip mined land is an 
Important issue In Kansas where the practice is in use and being expanded to new 
areas. 


GEOLOGIC PATTERNS ON ERT5 IMAGERY 

In satellite imagery of an orea such as Kansas with low relief, subtle 
geologic structure and extensive land-use, the geologic ground pattern that is 
most apparent is that ccused by stream patterns. Pattern and frequency of streams 
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ore important Indicators of the rock“type upon which they are developed- In 
addition, they are strongly influenced by geologic structure. Spatial frequency 
analysts of E RTS imagery lends itself to the study of drainage patterns by detecting 
those patterns which display a preferred orientation or spacing. Curves that result 
from such analysis may contain “signatures" that are attributable to basic geologic 
parameters. In addition, insight may be gained concerning structural trends by 
detecting stream orientations that may be influenced by joints and other structures. 
Such insight would be useful since little knowledge Is available concerning joint 
trends in Kcnsas and their significance . 

Use of spatial frequency analysis in the study of stream patterns should be 
guided by a knowledge of the manner in which stream patterns are expressed on ERTS 
imagery and how this expression vanes from place to place in the state os a result 
of changing rock-type, landuse, climate and natural vegetation. 

In Kansas, stream courses are generally expressed on ERTS imagery' in four 
different ways: 

1 . Riparian Vegetation 

In the eastern half of the state are many stream valleys that support denser 
end higher stands of vegetation than do the surrounding uplands. Cottonwoods, 
willows and other trees and shrubs which require large amounts of water thrive near 
streams that are perennial or contain water during much of the year. Such streams 
ore relatively easy to identify on ERTS- 1 images with the MSS5 band giving the best 
expression. On this band, riparian vegetation generally eppears much darker than 
surrounding fields and grasslands. In the western half of the state the increased 
dryness restricts this type of vegetation to only the major perennial streams such 
as the Arkansas River and the lower Smoky Hill. 

2 . Pure Topoqrcphic Enhancement 

Several areas in the state are largely uncultivated . For the most part, 
these are ranching areas that are covered with both natural and introduced grasses. 
Trees and shrubs are generally lacking even along streams in many of these areas. 

The absence of detractive patterns caused by fields and gross vegetation differences 
permit the enhancement of topography by differential illumination of slopes of 
varying orientation. This expression can be found on images covering the Kansas 
Flint Hills, as well as the Red Hills, Smoky Hills and dissected regions adjacent 
to larger streams in the High Plains. 
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Topographic Enhancement of drainage patterns has a much broader application 
In the winter at times of deep snow cover. The snow has the effect of giving the 
landscape more uniform reflecting properties by masking over areas of different 
crop and vegetation type. As a result, slope orientation becomes more critical 
in determining the amount of sunlight reflected to the ERTS-sensors. The lower 
sun angle in the winter serves to further enhance the topography. Thus an area 
in which stream patterns are not normally discernible may display them fully when 
snow-covered. 

3. Lond-Use 

Differing land-use between stream valleys and uplands can often accentate 
stream patterns. This can come about in a number of woys. One involves bottom 
land cultivation in an area of upland grazing and occurs in association with the 
major streams in the Flint Hills and other hilly areas, where the level fertile flood 
plains offer the most desirable farming areas. 

Another method by which land-use reflects stream patterns occurs in the 
western part of the state and is the opposite of the previously mentioned method. 

It is best displayed in the High Plains and dissected High Plains area where stream 
vaileys are often rough end lack flood plains. In this ccse the level upland areas 
offer the most ideal conditions for cultivation. The valleys which are usually too 
rough or rocky for farming are used as pasture. 

4. Direct Stream Expression 

In some situations, the actual stream beds can be delineated on ERTS imagery. 
This occurs in two ways. In the eastern portion of the state the larger streams 
usually contain water throughout the year and can often be discerned on the MSS 6 
and MSS 7 bands of ERTS imagery due to the low return of infrared energy from 
water bodies. Thus the lerger streams are often dark in comparison to the surrounding 
countrysides. 

In a different manner, cctual stream beds in the western and south-central 
part of the state are also expressed on ERTS imagery. In this situation, it is the dry 
stream beds which give o distinctive appearance. These dry streams are choked with 
sand, which highly reflects energy in the visible region and produces a bright 
appearance on MSS 4 and MSS 5 images which contrast well with the less bright 
appearance of surrounding fields and vegetation. 
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APPENDIX D. Fortran IV Programs 


This appendix contains a listing of the Fortran IV 
programs used in the analysis of the data obtained from 
the optical processing system. The programs are labelled 
and arranged in the following order; 

Program; 

1. DPSU DATA CALIBRATION PROGRAM 

Calibrates and normalizes the data obtained from the 
DPSU. Prints, punches, and plots the resultant data. 

2. DPSU DATA ANALYSIS --PARAMETERS /WEDGE CURVES 

Uses wedge data from program 1 to obtain parameters 
which describe features in the wedge data curve. Cal- 
culates these parameters for 4 sectors; wedges (1-8) , 
(9-16), (17-24), (25-32). Prints and punches resultant 

data . 

3. DPSU DATA ANALYSIS --SPATIAL FREQUENCY CURVES MODIFIED 
Calculates a modified spatial frequency curve by 
dividing each frequency curve obtained from program 

1 (point-by-point) by average frequency curve. Prints 
punches, and plots resultant data. 

4. DPSU DATA ANALYSIS— PARAMETERS/SPATIAL FREQUENCY CURVES 
Uses modified spatial frequency curve data obtained 
from program 3 to calculate parameters which describe 
features in frequency curve. Calculates these parameters 
for two frequency bands. Prints and punches the resultan 
data. 
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Ot 


C 

C 

C 

C 

c 

c 

c 

c 

c 


DPSl DATA CALIBRATION PROGRAM ( 05/25/73 ) 


PROGRAM 

INPUT — D(ItJ) (DATA FROM DPSU) “ " 

— N, (TOTAL NC. CF DATA SETS) 

N 1 , ( NO • OF FIRST DATA SET) 

OUTPUT-CALIBRATED RING, WEDGE 0 A T A- -PR I NT E 0 ,PU N C H E 0 V " 
PLOTTED 


LABELS - ■ 

L = NO * OF DATA SET 

A= RING AREAS ( 32 VALUES ) 

F= CALIBRATION FACTORS " ( " 3 2 ' ~"V A L U E $ )' - 
CORRECTION FOR WEDGES ) 


-TARE A /SENS ITT VI TY 


X ( 33 ) , 


DIMENSION 0 ( 70 , 64) ,6 ( 64 , 70 ) , CN (64) ,DL ( 64) , 01P (33 ) ,02P (33) , 

1A (32) ,F(32> ,GI (64) ,0T (64) ,KK ( 2 00 ) , LL (11 ) 

♦♦DIMENSIONS FOR O', B, SET' HERE FOR 70 OAT A” SETS** • 

DATA A/80,8, 69. 8, 202. ,344. , 5A8. ,707. ,10 5 0. ,1440. ,1920., 2510. ,3150. 
1,3910. ,4770. ,5730. , 69 50. ,8060. ,9390. ,10900. *12400., 14200., 16100. ,1 
100. ,20200. ,2 25 00., 25100. ,27500. , 30 50 0 7, 3 34 001 V 36 7 0 07,7 0 5 0 J7T4 3 3 0 
.,46400./, 

/ n i. 7 . nco * n>, . nr, 


28100 
30 
4F/ 
545 
625 
737/ 


00. ,20200. ,2 25 00 ., 25100. ,27500., 30 50 or, 33 4 0 0 .", 36 7 0 0 7,40 5 0 07743 3T 

,46400. /, 

1.052,1.070,1.043,1.058,1.061, 1.064,1.064, 1.054, 1.049, 1.0 
,1.037, 1.017, 1.008, 1. 008, 1.014, i. 0 1 7 , 1.008,170 06,1 . 0 0 0717 01 l,l70~ 

,1.031, 1.037, 1.034, 1.034, 1.040, 1.03^, 1.034, 1.0 37, 1.040, 1.037, 1.0 

An rq-i.nnTw.KH 


READ (5 , 10 0 ) N,Ni 
♦♦READ" OAT A USING FREE FORMAT ~ SUBROUTINE**" 
CALL FFMT (5,6) 



oocj o o o onn 


dg i j = i , e 4 

DO 1 1 = 1 , N 

1 0 (I, J> =9 ( J,I ) 

DC 40 K= 1,11 
40 Ll (K ) =K 

no 4 K=l, 33 

4 X (K) =K 

00 2 1 = 1 , N 

**RING DATA COM PU T A T I ON * * 

01 = 0 (1,1) /A (1) 

DO 3 J = 1 , 32 

ON 1 ( J) =0 (I ,J) / ( Ol * A ( J ) ) 
QT(J)=DN(J)*10.**10' 
OL(J)=ALOGlO(OT(J) ) / 1 0 • 

3 CONTINUE. 


**WE0GE DATA COMPUTATION** 

0 2 = 0 ( I * 3 3 ) *F C 1 ) 

DO 25 J = 3 3 * 64 

DI ( J) =0 (I , J) *F ( J-32> 

IF (01 (J) .GT. 02)02=01 ( J) 

29 CONTINUE ‘ “ 

00 50 J = 3 3 , 64 
50 ON (J) =01 { J ) /G? 


VALUES FOR PLOTS PUT INTO 01P , 02P 

OlP(33) = 0 .4 ' * : 

D2P(33> = 0 .4 
DC 7 J= 1 , 32 
01P(J >=QL (J) 



o o o 


7 02P(J)=GK(J+33) 
L = I + N 1 - 1 
K < ( L > - L 


**PUNCH* PRINT, PLOT, CALIBRATED DATA ; 

WRITE (4 3, 102) KK ( L ) , ( ON U > * J = 1 , 6 ) , LL ( 1 ) 

WRITE {43, 102) K K ( L ) , ( DN ( J ) , J = 7 ., 12 ) , L L < 2 ) _ 

WRITE (43, 102 ) K K ( L ) , ( CN t J ) , J = 1 3 , 1 8 ) , L L ( 3 ) 

WRITE (43, 102) KK(L),(CN(J),J-15,24),LL<4> 

WRITE (43, 102) KK(U»(DN(J),J = 25,3 0),LL<5)_ 

WPITE (43, 102) KK(L),(DM(J),J=31,36),LL(6) 

WRITE (4 3, 10 2) KK ( L ) , ( DN ( J ) , J = 3 7 , 4 2 ) , LL ( 7 ) 

WRITE (4 3, 102 ) KK (L ) , ( CN ( J) , J = 43 ,4a ) , LL ( 8 ) _ 

WRITE (43 , 102) KK (L>, <CN (J) ,J = 49,54> ,LL(9) 

WPITE (43,102) KK (L) , (CN (J) , J = 5 5 , 6 0 ) ,LL (10) 

WRITE (4 3, 112) KK(L> , ( DN ( J) , J=6l ,64) , LL (ID 

' WRITE (6,113) KK (L) 

WRITE (6, 114) 

WRITE (6,101) (!) (I, J),J= 1,64) ‘ 

' ' WPITE (6, 104) 

WRITE (6,105) 

WRITE (6,106) ((J,DN(J),DL(J),J,DN(J+32)),J=1,32) _ 

WRITE (6, 1 15) ' ’ ' ' ' 

WRITE (6,116) 

CALL PLT2 (2,X,01P,33,X,D2P,33) _ ' _ 

WRITE (6,107) " ‘ “ ■ 

WRITE (6*108) 

WRITE (6*110) ' ' 

WRITE (6,111) ' 

2 CONTINUE 

100FORMAT(1X,I4,I5) „ 

"101 FORMAT ( 6 X ,8 F 8 * 3 ) > 

102 FORMAT ( 1 X , 1 3 , 2X , 6E1 1 . 4 * 6X , I 2 ) 

104 FORMAT (1HQ,5X,49HR = NORMALIZED RING VALUE, W = NORM AL 1 7 ED WEDGE VALUE 



Cn 


1 ) 

105 FORMfl T (3X.-PING NO.',7X,*P',HX,-ALCG10CRK»/10.',6X,'HEDGF NC.*,6X 

It W ) 

10 6 FORMAT {6XtI2,5X,E11.4,5X,Fl!?.A,10X*I2»5X,F10.4) 

107 FORMA T ( 1 1 X, + - - + •-+-- + - - 1 \ + _ - _' 4 . i 

+ + — 1 • - 4- - - + — + - + — 1-- + — + • ) 

108 FORMAT (16 X, , 2',5X,'4*,5X,*5',4X,*8 , ,4X, , 10',4X, , 12 , ,4X, # i4*,3X, 4 16 
1 tAX, 18* 1 4 X * ' 2 0 ' * 3 X * * 2 2 ' , 4 X , * 2 4 ' # 4 X , ' 2 E * » A X * • 2 8 ‘ , 3 X , * 3 G ' V4X * ' 3 2“* )" 

lie • « 7 

111 FORMAT (• WEDGES (X)- 1 ,3X, *15 %'iOK', *35^1 2Xi' f 6' 0' , 7i'2x7 - ‘S5" r yr2*T 

1110 »liXt , 135* tllX, *160* t5X,* (OEGRE FS) " ) 

112 FORMAT(lX,I3,2X,4Ell.A t 28X,I?) 

113 FORMAT ( * 1 • t • OA T A SET NO. * » IX, 13 ) : 

114 FORMAT (“ RAW DATA-') 

• m < «J52i; U . MN ?.. R AN0 W ARE PUNCHED. COLUMNS ALQGIQ (RK)/10 . AND • 

iA ARE PLOTTED VS. LlNES/Ml. AND DEGREE ST"’ RE'S P"E C~tT\? tiLT* *1 -• 

116 FORMAT C WHERE K=iQ.**10.M LY ’ 1 

__ ST0P 



non n n, i nnonnnn 


c D P SU Hat A ANALYSIS -- PARAUETERS/WGDGE CURVES 

C 

C 

c PROGRAM 

C INPUT-- DNI(N), CALIBRATED, NORMALIZED* RING AND WEDGE DATA 

C OUTPUT -- W1-U3, W9-W 16* W17-W2 4, Wa5-W3J>* OVERALL 

C 1> AVERAGE CHAV) 

c 2) RMS DEVIATION FRO.) AVERAGE (R M$D> 

C 3) AREA _ 

c <*)ABoVE/HELOW AV E )--(+) AREA =? ( AREA) * <-> AREA* < ARAN > 

C 4) DATA POINTS ABOVE AVE) (PAKS) 

5) NUMBER OF 1 p E A K S * IM REGION (PEAKS ) _ 

(NOTE-- REGION 1= ( W3.-W8 ) * RED ION 2=(W9~W16)* REGION 3s < W 17- W24 ) , 
REGION 4 s(W25-H32) ) ... . _ 


1H> _ . _ ■ 

DIMENSION DM ( 6 4 , 80 ) * D IFF < 32 # 8Q),KK<80} UMV( 6Q),DAV1( 80)* 
1DAV2( 3 0 ) * DAy 3 ( 39>>DAV4( 80>iAREA< 00)»AKtAl< 80'),ARfcA2< «0)» 

.2 ARE A3 ( SCI >.* AREA4.C -SQ ) * ARAiJK 8 u ) A R A.U 2. ( SO 1 *_ ARAM 3 3 8 Q.i_*_A R AN4J 8 OJLl 

3P AK3 ( BOJrPAKK 3Q>*PAK2< 8f)),pAK3< 30)»PAK4< 80)»RMSD< 80). 
4RMSDK 80)»RMSO2( 3p ) * RMSD3 ( 8U > » RMSU4 ( 30) * A1 ( 8 0 } * PE AK 1 ( 8 n ) , 

5PEAK2(30)*PEAK3<30) .PEAK4<80). . 

DATA NO /QO/ 


DO 309 I = l.MO 

....... READ. DATA 

READ(S*100) (DN(N. I ) *N=1»60) 

100 FORMAT (6X# 6Eil * 4 ) 

RE AD (5 *113) CKK( l ) , <DN(N, I ) »N = 6i*64) ) 



113 FORMAT ( 1 X i 13 , 2 k, 41-11.4) 
c ■;}• •:» v «■ o o- # ■:> **, # 

C 

C DATA AVERAGE CALCULATION 

C ... 

SUM1=0.00 
SUM2=0 . 00 
SUM3=0. 00 
SUM4=0.QQ 
DO 21 N s 3 3 » 4 0 

. 21 .. SUM 1* SUM 1 . DMT N » 1 1 ___ 

DO 2? 14 = 41.48 

22 SUM 2 = SUM 2 ♦ DM < N » I 

DO 23 N = 49 ? 56. . . __ __ 

23 SUM3sSUM3 ♦ DM ( H , I ) 

DO 24 Ns57»64 

.24 SUM4 = SUM4 * D N ( N * I ) 

S U M = S U M 1 + SUM 2 + SUM 3 + 3UU4 
DAV(I) 3 SUM/32,0 
DAV1( I ) = S IJ M 1/8 , on 

D A V 2 ( I ) =.. SUM 2/6 , QQ..„ ; 

DAV3( I ) = SUM3/A. 00 
DAV4( I ) r SUM 4/9 . 0 0 

C 

C 

c data area/data peak calculation 

area ( I ) = o.oo 
areaki > = 0,00 

AREA2 (I) =.0.00 

A R E A 3 ( I ) = 0.0 0 
AREA4 ( I ) = 0.00 

: AR A Nl( I )._a. ’0..00 

ARAN2I I > = 0.00 
A R A N 3 ( I ) = 0.00 



A R A M 4 ( I ) = 0.00 
TMRS = L) A V C I ) 

P A K 1 ( T ) = 0.0 

PAK2< I ) = 0.0 .... 

P A K 3 Cl) = 0.0 
P AK 4 ( I ) = 0,0 

. PEAK 1 ( I) = 0.0 

P&AK2C I > = 0 . 0 
P E AK3 ( I) = 0.0 

.... PEAM(I) = 0.0 

C 

c 

SUMSQ = 0.00 . 

DO 30 N = 33,40 

DIFF = DN(N. I > - DA VI < I ) 

SQDIF = DIFF * DIFF 
SUMSQ = SUMSQ ♦ SQDIF 

IF(DN(N, I ) .GT.TURS) PAK1 (I) s PAKl(I) ♦ 1.0 

DIFF(N,I) = DN(M, I) - t)AV.< I 

IF(DIFF(N. I > ) 31,32,32 

31 ARANl (I) = ARAMl < I ) - UIFF(N,I) 

DIFFCM, I) r 0,00.. . 

32 AREAKI) = AREA1CI) + DIFFC|,I1 
Y = DQ(N, I.) 

Y2 = DMC.M-1-. I) , 

I F ( N . E Q » 3 3 ) Y2 = [)U<N + li I ) ' 

I F ( N , E Q . 6 4 ) GO TO 33 

- -Yls DN<N + 1.-1) ' t 

' . GO TO 37 

35 Yl = Y2 

- 37 - 1 F ( ( Y . GT * Yl ) . AUU . < Y . G T .. Y 2 )..) _ P. u A K 1 ( .1 .BEAJMjLU 

_ _ 30 CONTINUE 

SUMAV = SUMSQ / 8.00 
RMSDK I ) = SQRT(SUMAV) 




SUMSQ = 0.00 

DO 4 1 ) i | = 4 X , 4 d / 

D I FT = !)N(N. 1 > - ! 1 A V 2 ( I > 

SOD IF = DiFF * DiFF 
SUMSQ = SUMSQ + 3Q0IF 

IF( D\!(N # 1 ) .St. TMliS) P AK2 (I) = PaK 2(IJ’ ♦ 1.0 
nlFFCI, I ) = UM( il, I ) - !) A V ( 1 ) 

I F ( D I r f < N , I ) ) 42, 43,43 

42..ARAN2 ( I ) = ARA!J2(I.) I FF ( N f j )_ 

D IFF ( r| , I ) = n,0fj 

43 ARE A2 { I ) = AREA?(I> + UIFFC'MJ 

- - Y = DM ( N # I ) 

Y 2 = DU ( N-l » I ) 

I F ( N . E Q . 3 3 ) Y2 s DM ( M + l i I) 

lFtN.EQ.64) . GO TO - 45 _ . , . 

Y 1 = DM(U + 1, I ) 

GO TO 47 

. 45 -Yi = Y2 .. . . _ ] i 

47 IF ( (Y.GT.Yl).AND. ( Y . G T ,Y2>) PEAK2U) = P 8 A K 2 < I ) +. %.Q 
40 CONTINUE 

SUM A V = SUMSQ ./ -8 ,-00- — - - 

RMSD2U) = SQRTtSUMAV) 

C 

SUMSQ = 0,00 

DO 50 M s 49,56 

DIFF = DM (Mi 1 > - (JAV3U ) 

SQDIF s DIFF » DIFF 
SUMSQ = SUMSQ * SQDIF - 

IF(DN(N,I).GT.THRS> PAK3 (I) = PaK3(I) ♦ 1.0 

- . DIFF(N.I). = DN(!'J,I>. - DAV(I) 

IF(DIFF(N, I ) > 52,53,53 
52 A R A N 3 (I) = a R A M 3 ( I > - DIFF(N,I) 

D.IFFCN, I ) s .0,00 



53 A R E A 3 ( I ) = A f \ T: A 3 ( l J * n 1 f F t I » I ) 

Y = DUO* I ) 

Y2 = DN(ii-i» 1 ) .;._i 

I F t N . E Q « 3 3 ) Y2 = 11N(N+1# I) 
lF<N.FU,64> GO TU 55 

Yl = DM < N + l » 1) - - .... 

GO TO 57 
55 Yl = Y2 

57 IF((Y,GT.Y1>. AND . ( Y , GT , Y2.) ) P.E A K .3 ( 1 1_*_ P S A K 3J„L)_*_X.*..0„ 

50 CONTINUE * ; 

SUM AV = SUNSU / 8.00 7 

RMS D3 ( I ) = SQRT(SUMAV) 

SUMSQ = 0,00 
DO 60 N s .57,64 

. ... . DlFF = Q M ( N * I ) D AVAIL). . . .1 . 

SQDIF = niFF « a IFF 
SUMSQ = SUMS Q + SQDIF 

IF (DM (N, I ) .QT.THRS) -P.AK4 ( I ) = .... P A K 4 ( I I t Q _ 

D I F F ( N » I ) = DM ( N , I > - DAV< I ) 

IF (DIFF(N*I)) 62*63*63 

_62- ARAN4 (I) = ARAN4(I) - DlFFCN. f I.) . : — 

D I F F ( N * I) = 0,00 

63 aREA4(I) = AREA4(H +' DIFF( UIJ 

Y2 = DN ( N-l » I) ^ * 

I F ( N . EQ , 33 ) Y2 a DN<N*1, D 

1 F ( N. EQ . 64 ) _GO„ TO- .6 5 — — — 

Yl = DN(N + 1» I) 

GO TO 67 

. 65 Yl = Y2- -• — 

67 IF(<Y.GT..Y1>.ANU.(Y,GT.Y2>) PEAK4 ( I ) a RE AK4 ( I > + 1.0 
60 CONTINUE 

„ SUMAV =. SUMSQ / 8.00 

R M S D 4 ( 1 ) a SQRT(SUUAV) 



c 


aREA< 1> = AfUiAKI) ♦ A K £ A 'M I ) + ARliA3(I) + AKfcA4(I) 

PAKS(I) = PAkKI) + P A K 2 < I ) + PAK3CI) + P A K 4 ( X ) 

SUMSQ =0.00 

. on 39 N = 33,64 - 

SQDIF = D l F F { U , I ) * OIFFCJ# I ) 

SUMSQ = SUMS a + SQ 3 1 F 

39 CONTINUE ... __ ___ 

SUMAV = SUHSQ/32,0 
RMSD(I) = SORT (SUMAV ? 

3.09. CONTINUE .. - : 1 

C O « tMHHUHUHHUHHHUHJ ft 

C 

. WR I TF ( 6 1 101 ) . ... _ 

WRITE(6,l02> ( KK( I ) # DAV( I ) iDAVK I ) , DAV2{ I ) >UAV3< I ) , DAV4< 1 ) » 

11=1. M0) 

. . .. WR I T£ ( 6 . 103 ) . 

WRITF(6,1i) 4)(KK( I ) » ARE A ( I ) » AREA1( I ) # AR5A2 < I ) » ARE A 3 ( I >, ARfcA4{ I ) , 

II si i NO) 

W R I T E ( 6 » 1 0 5 ) - - - 

WRI TE( 6, 106) (KK( I ) . AGAMIC I ) , AR AN2 ( I > , AR AN3 < I ) » ARAN4 < H . I «1 * NO ) 
l-JRITE(6,107) 

• WRITIH(6, 1 0 3 ) C K K ( I ) , P AKS ( I ) * PAKKI ) iPAK2( I.U.PAK.3C U.RAKAtU., 

. 11 = 1. M0) 

WRITE(6,109) 

WRITE{6,110)(KK(I),RHSD( I > , RMS0.1 U ) , RHSD2 ( I ) , RMSD3 ( t) , RMSU< < I > , 

~ WRITE (6. 121) 

WRITE (6, 123) < KK < I U P E A K 1 { I ) , PE AK2 < I ) ,PEAK3< 1 ) , PEAK4 ( I ) , 

C • : 

WR I TE ( 43 . 4 01 ) ( KK ( I ) # QAV ( I ) , DAV1 ( I ) » DAV2 ( I ) > DA V3 < I ) , DA V4 <1 ) , 

_lRMSDU),RMSDl(I),RnSD2(I) l RMSD3<I)_,RMSn4(I)_t_ in^Q) 

401 FORMAT ( I 5 . 10F7 » 3 ) 

WR1TF( 43,402) < KK ( I) » ARE A ( I > , AKEAK I ) , AREA2 ( I ) # ARE A3 < I ) » ARE A4 < I ) , 
1.ARANK I ) , ARANZ< I ) * ARAN3C I ) » ARAN4C l ) » l=l#NO)_ 



n o o o o o 


o 

4* 


?i>X » 1 ( - ) ARE A 1» , 6X, M-> AREA 2 « p6x. •<- > AREA 3S6X, 
) 

» I3,nXfFl5.3*.Flf»,3»Fl5.3.Fl5t3J _ 


lax. ’PARS' . iuX , 1 P A K 1 * * 1 0 X , »PAK2 MOX, *PAK3 * * 10X # 
l3iFn . 3. » F 1 5 . 3 ». F* 1 5 , 3 i F.15_,.3i F15. .3 Jl 


402 FOR.UT< 15, 90.0) 

W R I T p ( 4 ,5 # 4 U 2 ) (Ki;(I).PAK5CI).HAKKI> t PAK2(I).PAK3U)*PAK4U># 

1 PEAKlt I) »PEAK2U > r.PEAh3f II »PEAK4 ( D . 1*1, NQ). 

101 FORMAT ( '1» .10X. 'AVC* *12X# 'AVElMlX, *AVE2'-*llX, » A VE3 * , 3 1 X , * A VE 4 • ) 

102 FORMAT ( »o*. I3^Fll,3»Fl5,3/Fl5 t 3»Flb#3#Fl f ).3) 

103 FORMAT ( »i » , loX, t a K H A TOT AL 1 » 5X i * ARE A 1»,9X»*AREA 2 »,9 X 
1. ’ AREA 3 ' , 9X . ' AREA 4 ' ) 

.104 FORMAT ( *o* » l3,Fn f 3,Ft5 f 3 f Fi5 ( 3fFl5.3,rx5. J) 

105 FORMAT ( » 1 ' 

1 ' ( - ) ARE A 4 

106 FORMAT ( * 0 

107 FORMAT ( ' 1 • 

I * PAK4 * ) 

.10-8 FORMAT { 'O' 

109 FORMAT ( '1 * 

I I RMS 4 1 ) 

110 FORMAT ( *Q» 

121 FORMAT ( * i • 

1 ' PE AK4 ' ) 

123 F 0 R M A T ( * o ' * 1 3 . U x — ,F 1 5 3 , F 1.5. ,.3 * F 1 5 .3 , F 1.5 . 3 ) 

» «. # u 

» «■ -tn>o » * # tt o fl. » >> # ^ » 


10X. 'RMS' «12X. * R, IS I’.iaX.'RMS 2>,10X,'RMS 3'.10X.' . 

I3iFll,3»Fl5,3»ri5,3» FIS.. 3 » F15 . 3.1 

25X. ' PEAK1 ' ,10X. 'PEAK2 ' » 10 X# ' PE AK3 * » 10 X. 


STOP 

END 
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C uHSL JATA Ail Al Y Si 3 - - SPATxAl. FREQUENCY uUR VES -MO D i F I E J 

o 

c 

0 

PROGRAM 

1) OBTAINS PL 1 = Dl/AVE ,, PL2 = ul/TH 
.2) PuOTS PlI 

lABElS 

G'i = CALIBRATED » / NORMALIZED DA T A FROM OPSU 
Trt = CALIBRATED , NORMALIZED", CORRECTED DATA FROM OPSU 
(NO iMAGE INPUT) 

AVE • = AVERAGE DATA SET OBTAINED FROM AlL SETS 
DN = CALIBRATED , NORMALIZED , CORRECTED , DATA FROM OPSU 

INPUT --..ON 1_ _ 

OUTPUT — PL1,PL2, PLI PLOTTED ’ 


OiM.ENSj.ON UN (64) ,Di (80 , 32) , Pul (3 3 ) , PL2 ( 32 ) , TH ( 32 > , 

1A VE ( 3 3 ) , SUM1 (32) , KK ( 1 48 ) , X ( 3 3) f LL ( 3 ) ,AVE0(33) 

DATA TH/ _____ i . QOE-.0 • 93.Er.2i 4,* 20E-3 * 3 .73E-*f , 

13.77E-5, 

23.tt2E-5,2.27E-5,l. 45E-5 ,1 .0iE-5,7.50E-6»7,Q6E-6 f l.Q4E-5,2.39E-5 

32 . 1 6 E - 5 , ,i • 7 0 E -5 , 1 • 31E-5 , 9 • 71 £ - 6 , 3 . D 4E-6 , 5 . 92 E-6 , 4 • 52 E-6 , 3 . 2 1 E- b 

42.29E-6,1.7lE-b,1.25E-6,a.a2E-7,b.28E-7,5.01E-7t4*GlE-7»3.45E-7 
52 .9 7E-7 ,2 .5 EE-7 , 2. 51 E- 7/ 

DATA NO/4/. . ._ 

INITI AlIZE SUM 1 (N) 

0 0301 N = 1,32 ' •" 

301 3UM1 IN) = fl • 0 0 
• .. 0 0910 1 =1, NO 



o O O i ■ '! OJU I j oi o o 


* * k E A U U A T A * * 

C 

REAc (5, 1U u ) (uN(N) t N = 1,6U> 

1 uO FORMA T ( 6 X ,6Ell . 4 ) 

REAu(5,li3) KK.(I) , (ON (N) ,N = 61, 64) 

113 F0kMATllX,i3,2X,4E11.4) 

** COMPUTE OAT A SET AVERAGE ** 

DO 201 N = 1,32 

01 ( I , N ) =. UN(N) 

201 SUMl(N) =■ SUMl(N) ♦ 01(1, N) 

910 continue 

uo . au 1 N = 1,32. 

AVE(N) = SUM 1 ( N ) /F L QA T ( NO ) 

A VE M = A VE (N> *10 ***10 
. A VE 0 ( N ) = AlOGIU IAVEM)_/1G. 0 _ . 

.301 CONTINUE 

A VEu (33) = A V E O (32) 

WRITE (6,504) 

...5 04 FORMAT Cl * , ■ DATA AVERAGE* ) 

WRITE (6,503) ( (N, A VEIN) ,AVE0(N) ) ,N=1, 32) 
503 FORMAT C *,6X,i3,9X,Ell.4,E11.4) 

** COMPUTE Pci, ,PL2 ** 

_ 00 901 I = 1 , NO 

DO 310 N - 1,32 
PlKN) - 01 ( I , N ) / A VE ( N ) 

P L 2 ( N ) =_01(1,N)/TH(N) 

310 CONTINUE 

Pul ( 33) = PL1 (32) 

* * P R I N T ,, P l’ot"T7 _ PUNCH ** _ 

WRITE(6,7ul) K K ( I )_ 



U A r A SET NO 


IX ,13) 


7 Cl FORMAT IT,' 

WRITE ( 6 » 5 u 1 ) 

50 X FORMAT ( •O^aX^N^lUX^Pwl^luX^PLa*) 

WRITE (6, 502M (N,P^1 (N> ,P^2 CN> > *N = 1,32) 

5 u2 FORMAT ( 1 *,6X,I3»oXtF7.3t5X,F7.2) 

OQ30GK=1,33 _ 

3 U Li X ( K ) = K ' 

UALu PlT2(1,X,Pu1,33,X,Pl1»33) 

WRITE ( b , 107) _ ; : 

"WRITE ( 6 , 1 G 8 ) 

DO 76 K = li3 
7b U(K) = K 

W RI T E ( 4 3 , 8ij J KK ( I) , (>Ll ( N ) ,N = 1 , 12) , L L ( 1) " “ " . 

WRITE (43, 80 ) KK(I),(PLl(Ni,N=13,24),LL(2) 

W RITE (4 3,82) KK 1 1 ) , ( PL 1 ( N i , H =25,32) , L L (3 ) 

8u FORMAT (IX ,13 , i X , 1 2 F 6 • 3 » 2 X , i 1 ) 

82 FORMAT (1X,I3»1X,8F6.3 ,2bX,Il) 

_9G1_ CONTINUE __ t . 

OAc.1 PLT 2 (1 , X , AVE j , 33 ,X , AVED , 33) 

WRITE (0,107 ) 

W RITE (6,108) 

10 7 FORMAT (1 1 X t * — + — + - - 1 - - + - - + - ♦ • -~+ -- i 4- -- + 1 - - + — +-- + — ' 

1 +-- 1 --+-+ — +-- + -- 1 -- + — + — +-+-* 1 --+--+“) 

_1Q8_F0RMAT (16X, *2* ,5X, *4* , 5 X , * b ■ , 4 X , * 8 * , 4 X , * 10 * , 4X . * 12 * , 4X , * 14 * * 3X , * 16_ 
1* ,4X , * 18* ,4X, *20 * , 3X, *22* ,4X,* 24* , 4 X , * 2 6 * , 4 X ♦ * 28 * ,3 X , * 30 * , 4X, * 32*’) 
STOP ' 

END ' 



n n n 


dpsu Data analysis 


parauetcrs/spaual frequency curves 


c 

c 

c 

c 

c 

c 

c 

c 

c 

C-- 

c 

c 

c . 

c 

c 

c 

c 

c 

c- ■ 


PROGRAM ' MOOIFI5D 

INPUT-- DaTA(!J) *DATa FRO | MODIFIED SPATIAL FREQUENCY CURVES 

— KK( I ) , UAT A SET M 1 , 

OUTPUT-- Rl-'*21.R22-R32. OVERALL 

1) average ( d a v ) 

2) RHS DEVIATION FROd AVERAGE... CRMSJJ) 1 

3) AREA 

<*) AUOVE/BELUW AVE,--(a)AREAs( AREA), (-)AREAs(AKAN) 

(**> AUOVE/3ELOH 1,0 — ( *> AREA* ( OARA - > ARE A= < DARA.N.) 

Ml DYNAMIC RANGE — (DYNR) 

3 ) SLOPE OF CURVE <R22-R32>--< SLOPE) 


(NOTE-- REGION 1 = <Rl-R21>* REGION 2 = (R22-R32)) 


DIMENSION DATA (32) ,KK( 32) , D AV< 32) # DAVl < 32 ) »DaV 2<32) , DIFF<32) » 
1ARANK 3 2 ) » ARAN2( 32) , AREA< 32 > , AR EAi ( 32 ) * AHC-AZ ( 32 ) , D I FF1 ( 32 > » 

2DARANK32.) . DaRAN2(.32) , D A R A 1 (.3 2 > * D A » A 2 .< .3 2 1 / » A R. A ( 3 2 ) DY.NRU.32-?.> 

3DYNR2(32) > D Y N R ( 3 2 ) * RMSDC 32 ) , RjlSDK 32 ) . RMSD* < 32) • SLQPE< 32 > • X < 11 > . 
4 Y £ 1 1 > . Al(32) 

DATA no/32/ ... . - .. - 

DO 900 I = 1 # MO 


** READ DATA „ 


RE AD £ 5 , 7 0 1 ) (DATA(M) » N = 1 » 24 ) 

701 F0RMAT(5Xi12F6-.3) _ .. ; 

RE AD (5, 703) KK( I ) , (DATA(N) .N*2S»32> 
703 FORMATUXt I3,lX,aFfi.3) 



uu c 


C DA ~ A A V' E£ ! v A 1 1 : w 

C 

S U M 1, = 0,00 
SUM2 =0.0 
SUM = 0,0 
aREai ( i > s n .on. 

AREA2 ( I ) = 0 . Q0 
ARAN1. C I > = 0,00 
ARAN? < I > a 0 . 00 
DARA1 ( I ) = 0.00 
DARA 2 C I ) a 0 . 00 

DARAN.KI) =' 0.0Q . . 

D A R A N 2 ( I ) =0.00 
DO 21 Ns 1 * 21 

21 sum i = sum ♦ data? n ) __ 

DO 22 N =22,32 

22 SUM2 = SUM2 + i) A T A ( N ) 

SUM = SIJM1 ♦ SUM 2 

_ . .. D A V ( I ) a SUM/32, Q A 

D A Vi ( I ) = SUN 1/? 1.0 • ■ 

DAV2< I ) = SUM 2/11,0 

* * R E G 1 CM i ( R 1 - R 2 1 > -- AREA ""6F 0 UR V § # R M $,7 dY F 

HU = D AT A ( 1 ) 

XL01 = DaTACD 
SUM SOI =0.00 

DO 31 N = l,21 ' 

, DIFTE = DATA(N) - DAVl ( I ) 

SQDIT s diffe * diffe 

SUMSQ1 = SUMSQ1. .+ . SQDIF 

DIFF(N) = DATA ( N) - DAV( I ) 

IFCDIFF(N) ) 32,33*33 

32_ A RANI < I ) = ARANlt I ) - DIFFJ.N.) ; 

DIFF(N) = 0.0 

33 AREAK I ) = AREAk I ) + U IFF C I) 



onn 


DIFFl(fj) = lJATAci) - 1.00 
IFCDITFK \) ) 35,36*36 
35 DARANKI) = DARANl(I) - DIFFl(N) 

D I F F 1 ( !J ) = P.0 ' . . 

■ 36 DARAKI) = DARAl(I) ♦ DIFFI(N) 

IF(DATACN) .GT.MI1) Hi 1 - HAT A <N * 

IF(D4TA(N).LT.XL01>XU01 = LMTAIN) 

31 CONTINUE 

DYNRK I ) = HU - XL01 

SUMAV = S U M S Q 1 / 2 i * 0 

RMSDKI) = SORT (SUMAV) » 

/ 

• ttRem-JN 2 CR22-K32>-~AR1 =a’ OF CURVEr RHS7lJYN"N*T 

H 1. 2 = PAT A ( 22 ).. - - - 

XL02 = D A T A ( 2 2 ) 

SUMSO? =0*00 

. .. DO 41 H = 22., 32 .... ' • — 

DIFFE = DATACJ) -DA-V2 ( I > 

SQDIF = DIFFE .» DIFFE 

SUHSQ2 = SUMSQ2 -+-.SQDIF.. * 

DIFT(M) = DATA(N) - DAV(I) 

I F ( D I F F ( N ) ) 4 2* 43*43 

42 ARAN2U) = ARAH2(I)-“ DIFF.FJ) _ 

DlFF(N) = 0,0 

43 AREA2U) = AKEA2U) ♦ D I F F { ' | ) 

DIFFHN) = D ATA ( p ) . 1.00 

IF(DIFFKM>> 45,46*46 

45 DARAN2U) = DARAN2U) - DIFFi(N) 

. . DlFFl(N) =0,00 - — 

46 DARA2CI) = DARA2(l> + DIFFl(N) 



I F ( DAT A ( [| ) . C,T . M T ;> ) :|I2 s 3 A T A < u ) 
IF<IU T A<W> .LT,XLU^>XL3^ = D A T A ( N ) 


M = N - 21 
Y ( M ) = QATA(N) 

41 CONTINUE 

DYNR2U) = HI? - X U 0 2 ’ 
SUMAV = SUMSuP /11.0 
RMSD2 ( I ) = SORT ( SIJM AV } 


DARA(I) = UARA1CI) + D ARA2 ( I > 

... AREA(I) = AREAidJ ♦ AREA3H.) ’ ! 

SUhSt) = 0.00 . . 

SUMSQ = SUMS Q 1 + SUMSQ2 

SUMAV = SUMSQ/32. 0 1 : 

RMSDd ) = SORT(SUMAV) 

_ HI = HU . _ _ 

If <HI2.GE.HIi) ill = HI 2 
XLO = XLOl 

I F ( XL02 , LE . XLOl ) . XLQ .. = _X L 02 , 

DYNR< I) = HI - XLO 

— DO SOI J s Dll - --- — ' 

801 X<J) “ FLOAT C .J ) / 1 0 * 0 

CALL SLOP<X.YilDC.SLOPE*COEF) , 

sloped) s . Slope ; , - 

900 CONTINUE 

„ WR I TE < 6 » 101 ) : __ ... 

WHITE! 6,102) (KKC I ) ,DAV( I ) * flAVl < I ) • D AV2 ( I ) * I *1 # NO ) 

WRITE <6, 103) 

_„WR I TE < 6, 1 tH >.( KK < I DARE A < I ) , AREAl < I 1 *k REA? J J ). ,_A R A N 1 .(JJ jlAKAN?. LLLi - 

II = I* NO) 

WRITE<6, 105) 

WRITE! o,lQ6) (KKC I ).# &AR A ( I ) , BAR Ai< I ) , DAR A2 (J ) * PARANK l ) » PA.HA.N2 (J. 



o o 


ii = 1,1m 
W R 1 T r ( 6 # 1 1 1 ) 

white Co. 112) ( KK ( I >,R4SU< I ? » RNSDl < I ) . RHSD 2 C I ) , 1 = 1 , NO.) 

W R 1 T 5 ( 6 . 1 0 7 ) 

WR ITE (o, lfld) ( KK < I ) » D YNR < I > .DYNRl< I > * PYNH2 * I > . I *1 • NO ) 

WR1TP( o,109) , 

write (6#no) < kk ( i > i slope ( i > , i = i.no> 

WRITE (43. 2 Jl) ( KK t I > .DAV( I ) .IJAVK I > .DAV2( I ) • SLOPE ( I ) . 1*1. NO) 

2 0 1 FORMAT ( I 5* 4F7 .3 ) ' 

WRITE< 43.203) <KK( I ) » A R E A < I) . AREAl ( I ) # AREA2 < I ) . ARAN1 ( I ) . AKAN2( I ) , 
1DARA( I ) , DARAK I ) ,l)ARA2( I ) .DARAWK-I ) . TJARAN2 < I ) » I *1 » NO ) 

203 FORMATS I5.10F7. 3) - - . - 

WRITE (4 3, 2 05) < KK ( I ) .rMSjjU ) ,RMSD1( I > ,RMS02< 1 ) - DYNRC 1 ) » DYNRl ( 1 ) , 

1 DY OR 2 ( I > , 1*1. NO) 

. 205 FORM AT ( I 5 , 6F 7 . 3.) 


101 FORMAT ( * 1 * ,10X, 1 AV^* *12X. »AVEl' #,UX. 1 AVE2’ ) 

~ .. 102 FORMATt »0 r » I3,Fn.3,Fl5.3,Fl5,i3) 

103 FORMAT ( * j, * , 1 OX * »AREA TOTAL 1 .5X. * AREA i'i’9X, * AREA 2»,9X. »<~>AREA 1» 
1.6X. * <-> AREA 2' ) 

- 104 FORMAT ( * q * > l3,Fli..3,ri5.3,Fl5 i 3*Fl3.3.Fl5.3)_ _ 

105 FORMAT ( * 1 1 #1QX» * i . 0 Area T GT Al » # 2X . ' AR E A 1*.8X*»AREA 2’,9x. 

1 * (-)AREAl* »6X. * (-) AREA 2') 

' 106 FORMAT ( 'O’ > 13. Fli* 3 . Fl5,3.Fl5i3rF.l5..3< F15 . 3) 

107 FORMAT ( ' i ' , luX , ’ 0 YOR V, ilX , ? ’jYrlRl * » 10X , •* DYNR2 ’ ) 

108 FORMAT ( ' 0 » . I 3 . Fu . 4 . F15 , 3 . Fj.5 , 3 ) 

109 FORMATf *i« ,10X. ’SLOPE OF CURVE- ( R1U-R32 1 ' l ' 

110 FORMAK ♦o’ .13. Flo. 5) 

111 FORMAT ( * 1 ' .lOX. 'RMS 1 * 12X* f R‘lS l’.iOX.'RMS 2’) 

112 .FORMAT ( ’0’ . 1 3 . Fll,,4 . Fl5_. 4 . F15 , 4 ) 


STOP 

END 



subroutine Slop (XiY.^u^l )pe,coef> 

DIMENSION X(i),Y<l> 

X*M = N 
XON'E = n . 

yone*o. 

XSQAr^O. 

XY = 0 . . . 

DO 3 K=i,N % 

XSOAR=XSQAH^X(K>*X<K> 

- X3NE=X0UE+X(K) _ 

YONEsYONE*Y(K) 

XY=XY+X<K)*Y(K> 

5 CONTINUE 

denom=xone*xone-xn*xS3ar 

SLOPES ( XOME*YONE-XY*XN)/DE-NOM 

- C = ( XY *XONfl~ YONE«XSQAR ) /DEMON 

ANUMsO . 

A D E N - 0 . 

— X A V = 9 . _ 

YAVsO . 

XN = N 

- -.DO 8 I =1 ,.N 

X AV = XAV + X ( I ) 

YAV = Y AV + Y ( I ) 

8 CONTINUE ___ __ 

XAV-XA V/XN 

YAV=YAV/XN 

DO . 9 I = 1 , ;j 

A NUM = AN } JM+ < Y ( I )-YAV)*<Y< I>-YAV) 
ADENSADEN + (X< I)-XAV)*(XU >"XAV) 

9-.C0NTINUE- 

COEFsSLOPE*SQRT( adEN/ANUM) ~ 

HRITEC6, 7) C,S'-OPK#COEF 

7 FORMAK3F10.3) . _ 

RETURN 

END 
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